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around  the  seal  ring  in  order  that  the  effect  of  waviness  remain  Indefin¬ 
itely.  A  new  device  for  creating  a  moving  wave  with  no  internal  moving 
parts  is  described.  Experimental  results  show  that  the  device  operates 
reliably,  rising  this  devicekJjesults  from  wavy  seal  experiments.,, show  that  a 
moving  wave  can  reduce  wear  by  at  least  a  factor  of  200  compared  to  a  flat 
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A  model  for  a  wavy  seal  is  presented  and  refined  using  the  Elrod  algo¬ 
rithm  for  determining  cavitation.  Predicted  results  for  the  wavy  seal  com¬ 
pare  favorably  to  experimental  results.  Optimun  wavy  seal  design  condi¬ 
tions  including  the  number  of  waves  are  established. 

Theoretical  and  experimental  results  are  given  on  the  effects  of 
radial  taper  and  high  temperature  operation  on  seal  performance.  Thermal 
rotation  is  evaluated  experimentally  and  predicted  theoretically.  A  con¬ 
verging  radial  taper  greatly  reduces  seal  friction  and  wear  during  initial 
operation.  Results  show  how  wear  progresses  across  the  face.  High  temper¬ 
ature  operation  reduces  seal  friction. 

A  new  concept  of  a  self-forming  offset-land  seal  is  described.  One 
ring  of  the  seal  continually  machines  an  optimum  shape  into  the  mating 
face.  Theoretical  and  experimental  results  are  given. 
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Mechanical  face  seals  are  used  in  numerous  applications  in 
Naval  machinery.  These  applications  range  from  propeller  shaft 
seals  to  boiler  feed  pump  seals.  In  such  equipment  the  mechani> 
cal  seal  plays  a  vital  role.  When  such  seals  fail,  repair  is 
costly  both  in  terms  of  lost  time  and  direct  costs,  so  any 
improvement  in  seal  life  and  reliability  would  be  of  significant 
benefit. 

As  more  advanced  equipment  is  designed,  it  is  sometimes  dif¬ 
ficult  to  achieve  desired  performance  in  more  severe  service 
environments  with  the  present  state  of  the  art  of  seal  design. 
Thus,  an  improvement  in  seal  technology  would  serve  this  impor¬ 
tant  application. 

The  immediate  objective  of  the  research  herein  is  to  further 
the  understanding  of  mechanical  face  seal  lubrication  phenomena. 
The  ultimate  objective  is  to  develop  the  capability  of  designing 
contacting  face  seals  having  a  longer  We,  greater  reliability, 
and  for  extreme  environments.  Thus,  the  objectives  of  this 
research  are  compatible  with  mechanical  face  seal  needs  for  Naval 
machinery. 
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Chapter  1 
Introduction 


Mechanical  Face  Seals 

In  applications  where  a  rotating  shaft  must  pass  from  one 
fluid  region  to  another,  contacting  mechanical  face  seals*  play 
the  essential  role  of  minimizing  the  transfer  of  fluid  between 
the  regions.  Applications  of  face  seals  range  from  water  pump 
seals  to  process  pump  seals  to  propeller  shaft  seals.  The  per¬ 
formance  and  reliability  of  contacting  mechanical  face  seals  are 
of  great  importance  for  any  type  of  equipment  where  minimal  leak¬ 
age,  high  reliability,  and  long  life  are  necessary.  Even  for 
equipment  where  these  factors  are  not  so  critical,  seal  failures 
and  short  seal  life  lead  to  high  operating  cost  due  to  down  time 
and  maintenance  cost. 

Even  though  mechanical  face  seal  technology  has  been  stead¬ 
ily  improving  over  the  past  several  decades,  further  improvement 
in  the  state  of  the  art  of  seal  design  would  be  most  beneficial. 
Although  seals  having  an  acceptable  life  and  reliability  can  be 
designed  for  many  applications,  further  improvement  in  seal  life 
and  reliability  would  result  in  significant  cost  savings  to  the 
user.  Also,  there  are  numerous  mechanical  face  seal  applications 
where  seal  loading,  reliability,  life,  and  leakage  requirements 
are  difficult  to  achieve  within  the  present  state  of  the  art. 
Examples  of  such  applications  are  seals  for  pumps  for  nuclear 
power  plants  and  seals  for  large  diameter  submarine  propeller 
shafts.  Additionally,  the  friction  losses  in  face  seals  repre¬ 
sent  a  significant  fraction  of  energy  consumed  for  pumping  pur¬ 
poses.  Within  the  present  state  of  the  art  of  seal  design,  it  is 
very  difficult  to  design  a  low  leakage  seal  that  also  has  a  low 
friction  loss. 

*The  class  of  low  leakage  face  seals  where  there  is  definite 
contact  and  wear  of  the  faces  as  opposed  to  some  hydrostatic  or 
hydrodynamic  designs  where  a  definite  clearance  is  maintained. 


The  main  barrier  to  the  advancement  of  the  state  of  the  art 
in  the  past  has  been  that  the  mechanics  of  seal  operation  were 
not  well  enough  understood  to  be  able  to  reasonably  anticipate 
seal  performance  as  a  function  of  design  parameters.  It  is 
essential  that  the  lubrication  mechanisms  between  the  faces  be 
well  understood.  In  this  work,  the  results  from  further  research 
into  hydrodynamic  and  hydrostatic  lubrication  mechanisms  in  face 
seals  are  reported.  Much  has  been  learned  about  these  mechan¬ 
isms.  Based  upon  both  theoretical  and  experimental  work,  there 
is  a  strong  promise  that  effective  use  of  hydrodynamic  and  hydro¬ 
static  mechanisms  can  now  be  made  to  dramatically  improve  con¬ 
tacting  face  seal  performance. 

Seal  Lubrication 

The  mechanical  face  seal  consists  basically  of  two  annular 
rings  which  rotate  relative  to  each  other  and  which  are  pressed 
together  by  spring  and  fluid  pressures  (see  Figure  1-1).  In  con¬ 
ventional  seals,  the  surfaces  that  rub  together  are  generally 
manufactured  as  flat  as  possible  initially  so  as  to  minimize 
leakage.  The  effective  gap  between  the  faces  is  ideally  quite 
small  (order  of  1  urn)  so  that  leakage  flow  across  the  faces  will 
be  quite  small.  The  difficulty  in  designing  a  mechanical  seal  is 
in  maintaining  the  gap  at  a  very  low  value  while  at  the  same  time 
providing  a  definite  lubricant  film  between  the  faces. 

The  load  that  must  be  supported  at  the  faces  of  a  mechanical 
seal  is  due  primarily  to  loading  caused  by  the  sealed  pressure. 
The  unit  face  load  can  be  expressed  as  some  fraction  B  of  the 
sealed  pressure  where  B  can  be  made  greater  or  less  than  unity  by 
geometry  selection.  The  load  support  at  the  faces  is  derived 
from  fluid  pressure  and  mechanical  pressure.  If  the  fluid  pres¬ 
sure  at  the  faces  is  large  enough  to  support  all  of  the  load, 
then  there  will  be  no  contact  and  no  adhesive  wear.*  If  none  of 

♦There  may  still  be  abrasive  or  corrosive  wear  even  if  the  sur¬ 
faces  do  not  touch. 
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Figure  1-1.  Mechanical  Face  Seal. 


the  load  is  supported  by  fluid  pressure,  the  load  must  be  carried 
by  mechanical  contact  and  the  wear  rate  will  be  large. 

In  practice,  conventional  seals  often  operate  at  one  of  two 
extremes.  At  one  extreme,  a  large  gap  will  be  created  by  hydro¬ 
static  or  hydrodynamic  pressure,  all  of  the  load  will  be  sup¬ 
ported  by  fluid  pressure,  and  the  seal  will  leak  a  large  amount. 
At  the  opposite  extreme,  the  gap  will  close  completely  and  leak¬ 
age  will  be  very  low.  Only  a  small  fraction  of  the  load  will  be 
carried  by  fluid  pressure,  and  wear  and  heat  generation  will 
increase. 

Based  on  the  above,  it  can  be  concluded  that  an  effective 
seal  should  operate  between  these  two  extremes— having  both  ade¬ 
quate  fluid  pressure  load  support  and  low  leakage.  The  seal 
should  operate  so  that  it  just  touches  to  minimize  leakage  but 
such  that  the  load  is  carried  by  fluid  pressure,  not  mechanical 
pressure.  To  do  this  requires  that  any  fluid  pressure  generation 
mechanism  used  to  provide  load  support  to  the  seal  must  be  very 
carefully  controlled.  At  present,  in  commercial  contacting  face 
seals  this  lubrication  is  left  primarily  to  chance,  and  often 
such  seals  operate  at  one  of  the  extremes  mentioned.  Quite  com¬ 
monly  such  seals  will  operate  in  the  low  leakage  condition  where 
wear  progresses  at  a  definite  rate.  Such  seals  are  quite  satis¬ 
factory  for  many  applications.  However,  when  pressure  or  speed 
requirements  are  increased,  then  this  mode  of  operation  may  lead 
to  rapid  failure  through  wear  or  heat  checking. 

In  this  research  program,  attention  has  been  directed  toward 
studying  the  effects  of  waviness  as  a  source  of  hydrodynamic  and 
hydrostatic  pressure  in  face  seals.  Particular  attention  has 
been  focused  on  hydrodynamic  effects  in  relation  to  roughness  and 
low  viscosity  applications  where  some  touching  is  expected.  Wav¬ 
iness  was  selected  among  the  various  mechanisms  available  because 
it  is  controllable  in  seals  and  a  better  understanding  could  lead 
to  improved  seal  designs. 


4 


Wavy  Face  Seal 

The  concept  of  waviness  is  that  the  film  thickness  varies  in 
some  fashion  circumferential ly  around  the  seal.  Generally  speak¬ 
ing,  film  thickness  may  vary  radially  as  well  as  tangentially. 

h  =  h(r,0)  .  (1-1) 

In  the  present  work  interest  is  focused  upon  film  thickness 
shapes  of  the  following  functional  form 

h  =  hQ  +  f1(r)  cos  ne  +  f2(r)  sin  n6  .  (1-2) 

At  any  particular  radius  r  the  film  shape  is  periodic  with  n 
waves  around  the  seal  and  is  therefore  wavy.  However,  film  shape 
can  also  vary  in  some  general  manner  with  r. 

If  f^(r)  a  const  and  f2(r)  =  const,  then  the  faces  are 
always  radially  parallel.  This  component  of  film  thickness  var¬ 
iation  is  commonly  termed  waviness.  If  the  f^(r)  ^  const,  then 
the  faces  are  not  in  general  radially  parallel.  The  f . ( r)  are 
referred  to  as  tilt.  Thus,  the  film  thickness  shapes  of  interest 
are  combinations  of  waviness  and  tilt.  Since  at  any  radius  the 
film  thickness  is  wavy,  the  combination  of  waviness  and  tilt 
defined  above  will  also  be  called  a  wavy  film  shape. 

The  reason  for  choosing  film  shapes  as  described  by  equation 
(1-2)  as  a  subject  for  study  is  that  these  shapes  can  conven¬ 
iently  be  generated  by  planned  mechanical  distortions  in  a  seal 
ring,  and  the  shapes  also  include  common  modes  of  unplanned  dis¬ 
tortion  found  in  operating  seals.  For  example,  generally  seals 
undergo  a  uniform  tilt  due  to  pressure  and  thermal  deformation. 
When  rings  are  loaded  by  any  non  ax i symmetric  load  they  become 
wavy  as  described  by  equation  (1-2). 

Now,  for  the  sake  of  illustration,  assume  a  seal  has  a  wavy 
film  thickness  shape  given  by  equation  (1-2)  where  n  =  2.  Assume 
a  ring  is  pressed  against  a  perfectly  flat,  moving,  mating  ring. 
For  n  =  2  there  will  generally  be  2  or  4  regions  where  film 
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thickness  is  smallest  (only  2  are  shown)  and  where  contact  can  be 
expected  to  occur  (see  Figures  1-2  and  1-3).  In  general,  a  fluid 
pressure  distribution  develops  according  to  the  Reynolds  equa¬ 
tion.  Both  hydrostatic  and  hydrodynamic  components  are  included. 
Angular  rotation  of  the  mating  face  causes  pressure  to  build  up 
in  the  converging  regions  as  the  minimum  film  thickness  region  is 
approached.  In  diverging  regions,  cavitation  will  occur  as  the 
pressure  attempts  to  become  negative.  For  an  outside  pressurized 
seal  the  cavity  may  appear  as  shown  in  Figure  1-2  depending  upon 
the  f -,* ( r ) .  The  fluid  will  flow  in  streamers  across  the  cavity 
and  a  full  film  will  again  develop  at  the  downstream  cavity 
boundary. 

The  problem  to  be  solved  is  to  determine  the  fluid -and 
mechanical  pressure  distribution  for  the  given  configuration. 
Given  the  pressure  distribution,  then  load  support  and  leakage 
can  be  calculated,  and  the  effect  of  various  parameters  can  be 
studied.  Given  this  tool,  then  waviness  shapes  which  provide 
improved  performance  can  be  found. 

During  the  first  portion  of  this  research  program  only  ra¬ 
dially  parallel  face  waviness  was  considered.  Much  effort  was 
spent  in  developing  solution  methods  and  studying  predicted 
results. 

As  background,  the  pressure  distribution  for  wavy  parallel 
faces  has  been  solved  by  various  methods  by  Findlay  [1],  Pape 
[2],  and  Stanghan-Batch  and  Iny  [3]  for  perfectly  smooth  faces. 
These  results  show  that  even  a  small  waviness  produces  sufficient 
load  support  for  complete  liftoff.  However,  heavily  loaded  or 
low  viscosity  seals  do  not  operate  with  full  fluid  pressure  load 
support  and  complete  separation.  A  definite  wear  results.  Sur¬ 
faces  of  such  seals  contact  during  operation,  and  pressure  dis¬ 
tribution  is  affected  by  interactions  with  surface  roughness  (see 
Figure  1-3).  Seals  of  this  type  operate  in  a  mixed  friction 
regime.  Fluid  pressure  load  support  due  to  waviness  may  provide 
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a  significant  fraction  of  the  load  support,  but  asperity  contact 
must  provide  the  balance. 

As  a  starting  point  for  this  research  program,  these  impor¬ 
tant  effects  were  included  in  a  wavy  seal  lubrication  model.  In 
the  first  annual  report  for  this  project,  reference  [4],  this 
more  general  problem  was  solved  using  a  one-dimensional  theory. 

In  the  second  annual  report,  [5],  the  much  more  complex  two- 
dimensional  solution  to  the  above  problem  was  solved.  The 
effects  of  waviness,  roughness,  asperity  contact,  wear,  cavita¬ 
tion,  and  elastic  deflection  were  included  in  this  model.  Using 
this  model,  predictions  were  made  for  the  relative  wear  rate, 
friction,  and  leakage  as  a  function  of  roughness,  waviness, 
speed,  size,  pressure,  viscosity,  and  material. 

A  number  of  conclusions  were  reached  based  on  these  first 
two  annual  reports: 

1)  The  effects  of  roughness  on  hydrodynamic  lubrication  are 
not  completely  understood.  Certain  fundamental  ques¬ 
tions  remain  concerning  the  roughness  model  used. 

2)  As  to  the  potential  of  utilizing  hydrodynamic  effects 
caused  by  parallel  face  waviness  to  advantage  by  design, 
the  results  show  that  wear  rate  and  friction  can  be 
greatly  reduced  while  maintaining  leakage  at  acceptable 
levels. 

3)  While  a  comparison  of  predicted  results  to  experimental 
results  given  in  the  literature  is  generally  good,  data 
contained  in  the  literature  is  incomplete,  so  more  com¬ 
plete  experimental  data  are  needed  for  comparison. 

4)  In  low  viscosity  or  heavily  loaded  applications  where 
some  touching  is  expected  to  occur,  waviness  will  wear 
away  with  time  and  any  benefit  derived  will  be  lost. 

5)  Based  upon  data  for  some  commercial  seals  and  using  the 
model,  it  was  determined  that  there  was  insufficient 
accidentally  caused  waviness  to  produce  significant 
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hydrodynamic  effects  in  water.  One  cannot  generalize  to 
say  that  such  effects  do  not  occur  in  commerical  seals. 
However,  the  model  is  provided  so  that  the  question  can 
be  answered  on  a  case  by  case  basis. 

Item  1)  was  treated  extensively  in  the  third  annual  report 

[6] .  Even  after  this  analysis  certain  fundamental  questions 
remain  concerning  how  to  deal  with  roughness  in  lubrication  prob¬ 
lems.  A  contribution  to  the  literature  on  the  subject  was  made 

[7] .  This  thorough  analysis  led  to  conclusions  allowing  certain 
simplifying  assumptions  discussed  in  the  previous  annual  report 

[8] . 

Item  2)  was  also  treated  extensively  in  the  third  annual  re¬ 
port  [6].  A  methodology  for  the  design  of  a  wavy  face  seal  was 
developed  and  applied.  Theoretical  results  showed  large  reduc¬ 
tion  in  friction  and  wear  rate  compared  to  conventional  designs 
whereas  leakage  could  be  controlled. 

Concerning  Item  3),  the  second  and  third  annual  reports 
[5,6]  describe  a  test  apparatus  designed  to  test  the  wavy  seal 
theory.  This  apparatus  has  been  in  operation  for  more  than  two 
years  and  many  tests  have  been  conducted.  These  test  results  are 
reported  in  Chapter  2  and  compared  to  theory  in  Chapter  4  of  this 
report.  The  previous  annual  report  includes  additional  results 
[8]. 

Early  in  the  test  program  it  was  observed  that  the  type  of 
waviness  which  can  be  practically  applied  is  not  of  the  radially 
parallel  type.  Waviness  generally  consists  of  alternating  tilt 
plus  radially  parallel  waviness  as  described  by  equation  (1-2). 
8ased  on  these  considerations,  a  new  model  for  predicting  per¬ 
formance  was  developed  and  appears  in  the  previous  annual  report 

[8].  In  the  present  report,  significant  improvements  on  the 
model  have  been  made,  and  a  comparison  to  experimental  results 
has  been  made. 
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Concerning  Item  4)  above,  a  solution  to  this  problem  was 
first  proposed  in  the  first  annual  report  [4].  It  was  proposed 
to  move  the  waviness  slowly  around  the  seal  so  that  whatever  wear 
occurred  would  be  uniformly  distributed.  Then  the  shape  of  the 
wave  would  be  preserved  and  tests  using  a  constant  wave  could  be 
made.  The  concept  is  illustrated  in  Figure  1-4.  This  concept 
was  incorporated  into  the  test  apparatus  and  is  described  in 
detail  in  references  [5]  and  [6].  In  this  report,  a  new  concept 
to  move  the  waviness  with  no  internal  moving  parts  is  described 
in  detail. 

In  summary,  in  this  report  additional  experimental  results 
using  waviness  are  presented,  the  wavy  seal  model  is  further 
improved,  and  theory  and  experiment  are  compared.  The  results 
using  a  new  concept,  that  of  a  self -generating  seal  profile,  are 
reported.  Results  from  tests  of  the  effects  of  radial  taper  and 
high  temperature  environment  are  also  reported  and  compared  to 
theory. 
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Distorted  Ring  Under  Pressure 
After  Wave  Has  Moved 


Figure  1-4.  Moving  Waviness  Concept. 


Chapter  2 

Experimental  Program  and  Results 


Since  the  time  the  previous  report  [8]  was  written,  several 
series  of  tests  have  been  conducted  using  the  seal  test  appara¬ 
tus,  several  improvements  have  been  made,  and  a  new  type  of  wavi- 
ness  device  has  been  designed,  fabricated  and  tested.  These 
activities  will  now  be  described  in  detail. 

Test  Apparatus  Modification 

During  the  past  year  several  minor  improvements  were  made 
which  are  not  related  directly  to  a  test.  First,  an  electronic 
tachometer  was  installed.  A  gear  was  mounted  on  the  drive  shaft 
adjacent  to  a  proximity  probe  such  that  the  output  voltage  under¬ 
goes  a  cycle  as  each  tooth  passes.  This  signal  is  fed  into  a 
pulse  counter  card  in  the  HP  6940  multi  programmer.  Using  a 
counting  period  of  400  ms,  controlled  by  a  timer  card,  a  subpro¬ 
gram  in  the  HP  9835  calculates  the  rpm  directly  and  displays  the 
information  during  a  test. 

Second,  a  reversing  control  was  installed  on  the  main  drive 
motor.  Whereas  previously  a  torque  zero  was  found  by  hand  crank¬ 
ing  the  system  first  forward  and  then  backward,  this  procedure  is 
now  carried  out  using  the  drive  motor  itself  to  provide  forward 
and  backward  rotation  at  low  speed.  This  procedure  has  the 
advantage  that  it  is  quicker  and  more  consistent  than  the  hand 
method . 

Third,  a  fixture  was  designed  to  provide  a  support  for  the 
pressure  vessel  when  the  vessel  is  installed.  Formerly,  instal¬ 
lation  of  the  vessel  was  somewhat  awkward  and  required  two  peo¬ 
ple.  The  device  is  designed  to  be  removed  from  contact  with  the 
vessel  after  installation  so  that  torque  readings  are  not  affect¬ 
ed.  Although  fabricated,  the  device  has  not  been  installed  at 
this  time. 
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Waviness  Drive 

As  discussed  previously  [8],  moving  the  wave  by  physically 
rotating  the  pressure  pockets  relative  to  the  seal  ring  has  the 
disadvantages  that  the  0-rings  slide  tangentially  and  are  subject 
to  wear  and  that  a  turning  mechanism  is  required  within  the  seal 
assembly.  As  will  be  shown  in  Chapter  3,  the  sliding  and  the 
turning  mechanism  can  be  eliminated  by  providing  three  waviness 
causing  pressures,  each  varying  sinusoidally  with  time,  to  three 
separate  sets  of  pressure  pockets  under  the  seal.  By  superposi¬ 
tion  of  deflections,  this  causes  the  wave  to  move. 

A  device  to  provide  these  pressures  has  been  designed,  fab¬ 
ricated,  and  tested  on  a  three-wave  seal.  Figure  2-1  shows  the 
waviness  drive  assembly.  The  three  cams  @  are  driven  by  a  low- 
speed  drive  (the  same  drive  as  used  previously  for  the  waviness 
drive)  at  5.3  revolutions  per  hour.  Each  cam  is  spaced  120°  from 
the  other  two.  The  cams  are  circular  in  shape  and  move  the  flat 
face  followers  ®  .  The  circular  shape  of  the  cam  displaces  the 
follower  in  a  simple  harmonic  motion.  Then  the  linear  spring  © 
produces  a  simple  harmonic  sinusoidal  force  as  a  function  of 
time.  These  forces  act  on  pistons  ®.  The  pistons  are  sealed' 
off  by  0-rings. 

At  the  bottom  of  each  piston  cylinder  is  a  valve  assembly 
shown  in  detail  in  Figure  2-2.  Oil  hydraulic  pressure  is  sup¬ 
plied  at  point  ©  .  If  the  piston  is  too  low  (meaning  that  the 
pressure  in  chamber  ©  is  lower  than  it  should  be  as  dictated  by 
the  spring  force),  then  the  needle  valve  seat  at  ©  opens,  ad¬ 
mitting  high  pressure  oil  into  the  chamber  ©.  If  the  pressure 
is  too  high,  then  the  piston  will  raise  and  open  the  needle  valve 
at  ®,  thus  releasing  pressure  from  the  chamber.  In  this  manner 
the  pressure  in  the  chamber  is  controlled  to  follow  the  spring 
force  so  that  the  chamber  pressure  (connected  to  the  seal  pres¬ 
sure  pockets)  also  varies  sinusoidally  with  time. 
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Section  on  B-B 


The  mechanism  in  essence  consists  of  three  programmed  pres¬ 
sure  regulators.  One  major  difference  between  this  device  and  a 
normal  pressure  regulator,  however,  is  that  both  an  increasing 
pressure  and  a  decreasing  pressure  can  be  provided  as  needed  as 
the  required  pressure  both  increases  and  then  decreases  with 
time. 

As  first  conceived  and  fabricated  the  valve  assembly  in  Fig¬ 
ure  2-2  was  not  installed.  Instead,  each  of  the  three  chambers 
0  and  the  lines  to  the  seal  were  filled  with  grease.  Theoreti¬ 
cally,  this  will  provide  the  same  result,  given  an  incompressible 
system.  However,  it  was  found  that  this  approach  did  not  work 
because  of  uneven  thermal  expansion  of  various  parts  of  the  sys¬ 
tem  and  because  of  deflection  interactions  at  the  seal  itself. 
That  is,  as  one  pressure  increased,  the  seal  would  deflect  and 
this  deflection  would  increase  the  volume  of  the  adjacent  pres¬ 
sure  system.  The  system  was  unable  to  compensate  for  these 
errors  so  an  active  system  using  the  valves  was  installed. 

A  flow  diagram  of  the  waviness  drive  as  connected  to  the 
seal  is  shown  in  Figure  2-3.  Hydraulic  pressure  is  maintained  by 
an  air  driven  oil  pump.  The  accumulator  reduces  pulsation  caused 
by  the  pump.  The  oil  is  supplied  to  the  valves  as  described. 

The  three  regulated  pressures  are  transmitted  to  three  sets  of 
pockets  under  the  seal  as  shown.  Leakage  flow  through  the  valves 
is  returned  to  the  sump. 

The  operation  of  the  waviness  drive  was  verified  experimen¬ 
tally  by  measuring  pressure  as  a  function  of  time.  Figure  2-4 
shows  how  the  system  responded  using  high  pressure  nitrogen  as  a 
supply  rather  than  oil.  The  curve  has  a  sinusoidal  shape  with 
time.  The  flattening  at  the  top  and  bottom  is  caused  by  friction 
in  the  system  and  the  error  introduced  by  the  valve  opening 
force. 

Figure  2-5  shows  the  time  response  using  high  pressure  oil. 
Although  the  general  shape  is  clearly  sinusoidal,  there  are 


Figure  2-3.  Flow  Diagram  for  Waviness  Drive  System. 


Pressure  Measured  Without  Carbon 
Ring  in  System  and  Using  N.,  lias 
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Figure  2-4.  Waviness  Pressure  with  Time  -  Gas  Supply. 


Pressure  Measured  with  Carbon 
Ring  in  System  and  Using  1)0  Uil 
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Figure  2-5.  Waviness  Pressure  with  Time  -  Oil  Supply. 


clearly  jumps  in  the  pressure.  Again,  these  are  believed  to  be 
caused  by  0-ring  friction  coupled  with  the  higher  viscosity  of 
the  oil.  It  is  expected  that  some  refinements  will  be  made  to 
the  system  to  reduce  these  errors,  although  it  Is  thought  that 
there  is  little  effect  on  actual  seal  operation. 

Even  though  the  time-dependent  behavior  of  the  pressure  is 
better  using  gas  than  oil.  It  was  decided  to  use  oil  for  the  sys¬ 
tem  because  an  oil  supply  system  is  simple  and  there  are  no  prob¬ 
lems  with  corrosion.  Long-term  reliability  is  expected  to  be 
better . 

Using  this  device  one  500-hour  test  has  been  completed  (Test 
110  in  Appendix).  No  difficulties  were  encounted  in  the  test; 
the  device  operated  reliably.  A  comparison  of  test  results  for 
the  seal  shows  that  seal  performance  is  quite  close  to  that  ob¬ 
tained  using  the  previous  moving  wav i ness  device.  Thus,  it  has 
been  demonstrated  that  the  waviness  drive  concept  works.  Given 
the  inherent  reliability  of  the  waviness  drive  system,  it  is  ex¬ 
pected  that  it  will  operate  indefinitely,  whereas  the  sliding  0- 
ring  approach  would  be  expected  to  fail  because  of  0-ring  wear. 
(The  500-hour  test  using  sliding  0-rings  did  show  signs  of  deter¬ 
ioration.) 

Friction  and  Wear  Test  Apparatus 

In  comparing  the  results  of  certain  experiments  to  the  pre¬ 
diction  of  theoretical  models,  it  has  become  clear  that  the  mech¬ 
anical  coefficient  of  friction  between  carbon  and  tungsten  car¬ 
bide  is  not  a  constant  but  varies  widely  depending  on  contact 
pressure,  temperature,  and  the  amount  of  water  present.  Similar 
inconsistencies  are  observed  for  wear  in  that  the  wear  coeffi¬ 
cient  used  in  the  simple  wear  law  is  apparently  not  a  constant. 

What  is  needed  is  a  set  of  fundamental  experimental  data  re¬ 
lating  friction  and  wear  to  various  operating  conditions.  A 
search  of  the  literature  has  shown  that  such  data  does  not  exist 
for  tungsten  carbide  versus  carbon  operating  in  water.  It  is 
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known  that  mechanical  friction  coefficient  is  sensitive  to  the 
presence  of  water,  so  dry  friction  measurements  are  not  very  use¬ 
ful.  Such  data  does  exist  for  these  materials  in  a  mechanical 
seal  configuration.  However,  this  data  is  of  little  use  because, 
as  is  well  known  now,  there  is  no  way  to  separate  the  actual 
mechanical  contact  pressure  from  the  fluid  pressure  in  a  seal,  so 
any  estimate  of  mechanical  friction  coefficient  based  on  seal 
operation  is  a  guess  at  best.  Mechanical  friction  coefficient 
must  be  measured  outside  of  the  context  of  a  seal  where  fluid 
pressure  load  support  can  be  eliminated  or  at  least  limited  to 
microasperity  lubrication  caused  by  surface  roughness. 

Given  that  such  data  does  not  exist  and  is  important  for  the 
development  of  refined  seal  models,  the  design  of  a  friction  and 
wear  test  apparatus  has  been  undertaken.  Figure  2-6  shows  the 
design.  The  friction  pair  consists  of  three  carbon  cylindrical 
plugs  @  being  pressed  against  a  rotating  W-C  ring  ®  .  Three 
plugs  are  used  to  provide  symmetry  of  load  which  simplifies  the 
design.  The  plugs  are  small  enough  such  that  hydrodynamic  and 
hydrostatic  effects  should  be  minimized  relative  to  the  applied 
mechanical  load.  Load  is  applied  by  air  cylinder  ©  acting  on 
the  end  of  shaft  ©  and  in  turn  onto  the  mounting  ring  for  the 
three  plugs.  Shaft  ©  is  designed  to  easily  bend  to  compensate 
for  uneven  wear  of  the  plugs  and  wobble  of  the  face  of  the  rota¬ 
ting  ring.  Shaft  ©  transmits  the  friction  torque  to  reaction 
arm  ©  on  which  strain  gages  are  mounted  to  indicate  torque. 

The  entire  assembly  is  housed  in  a  pressure  vessel  ©  so  that 
the  contacts  can  be  flooded  with  water  and  the  water  can  be 
heated  under  pressure  up  to  400*F.  The  pressure  is  contained  at 
the  drive  shaft  end  by  mechanical  seal  ©  and  at  the  torque 
shaft  end  by  a  small  diameter  0-ring.  The  0-ring  is  calculated 
to  introduce  only  a  small  error  into  the  torque  readings.  A  var¬ 
iable  speed  drive  motor  is  connected  to  the  drive  shaft. 

At  the  present  time  detail  drawings  of  the  apparatus  have 
been  completed,  and  fabrication  will  begin  soon.  The  first 


series  of  tests  will  be  conducted  at  room  temperature  under  both 
dry  and  flooded  conditions  at  various  speeds  and  mechanical  con¬ 
tact  pressures.  Both  friction  and  wear  will  be  measured. 

Radial  Taper  Tests 

A  series  of  tests  was  conducted  for  the  purpose  of  evaluat¬ 
ing  the  effects  of  various  initial  radial  tapers  on  seal  per¬ 
formance,  both  short  term  and  long  term,  and  for  the  purpose  of 
comparing  to  results  based  on  a  theoretical  model  developed  pre¬ 
viously  in  relation  to  this  project.  The  model  considers  the 
effect  of  initial  taper,  thermally-caused  taper,  and  wear  on  seal 
performance  [9]. 

As  series  of  14  tests  beyond  those  reported  previously  [8] 
were  conducted.  These  tests  were  conducted  for  a  balance  ratio 
B  =  0.75  so  the  original  seal  design  was  modified  as  shown  in 
Figure  2-7.  The  purpose  of  using  B  *  0.75  was  to  be  able  to 
achieve  complete  liftoff  at  high  convergent  tapers. 

Snapp  and  Sasdelli  [10]  produced  a  radial  taper  by  lapping 
seal  rings  on  a  commercial  lapping  machine  and  manipulating  the 
machine's  operating  conditions.  Experience  shows  that  this  meth¬ 
od  is  effective  but  can  be  slow  and  unpredictable.  The  needed 
manipulation  is  that  of  moving  the  lapping  rings  in  toward  the 
center  of  the  lapping  plate  or  out.  Moving  the  lapping  rings 
inward  causes  the  lapping  plate  to  wear  convex,  thereby  producing 
a  carbon  ring  with  a  converging  radial  taper.  Moving  the  rings 
outward  has  the  opposite  effect.  This  method  takes  time  to  wear 
the  lapping  plate  to  the  desired  shape  and  it  is  not  easily  con¬ 
trolled. 

The  two  devices  shown  in  Figure  2-8  have  proven  to  be  both 
reliable  and  easy  to  use  for  producing  a  uniform  moment  about  the 
ring  centroid  causing  the  faces  to  either  converge  or  diverge. 
Device  A  is  used  to  create  a  net  divergent  taper  by  placing  it 
into  the  seal  ring  close  to  the  sealing  face.  It  is  then 
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pressurized  and  the  resultant  profile  becomes  convergent.  The 
seal  is  then  lapped  flat,  the  device  depressurized  and  the  result 
is  a  divergent  taper.  Device  B,  an  expanding  mandrel,  is  placed 
at  the  extreme  opposite  end  from  the  face  of  the  seal.  Because 

of  the  slots  in  the  device,  when  the  bolt  is  tightened,  the 

wedges  move  inward  expanding  the  apparatus.  The  0-ring  is  used 

to  help  even  the  load  around  the  circumference.  This  expansion 

causes  a  divergent  taper  which  when  lapped  flat  and  the  device 
removed  results  in  a  net  converging  taper. 

The  test  procedure  used  for  the  tapered  seal  tests  is  simi¬ 
lar  to  that  reported  previously  [8].  The  Appendix  (Tests  64-68, 
75,  76,  97-101)  contains  performance  curves  for  the  experiments 
conducted.  Table  2-1  is  a  summary  of  results.  Test  64  was  run 
with  a  carbon  seal  which  had  an  initial  slope  of  1013  um/m  con¬ 
vergent  taper.  It  was  expected  that  torque  would  be  quite  low 
and  leakage  high.  Table  2-1  shows  that  the  torque  was  0.05  N*m 
and  the  leakage  very  high  at  58  cm3 /min. 

Test  64  in  the  Appendix  shows  an  erratic  operation  with  a 
sudden  increase  in  torque  and  seal  face  temperature  approximately 
4.5  hours  after  start  and  a  sudden  decrease  at  approximately  25 
hours  after  start.  This  closely  corresponds  to  18:45  and  05:30 
clock  time.  These  jumps  are  also  seen  in  Test  65.  The  periodic 
nature  exhibited  called  for  a  closer  examination.  Several  tests 
were  made  with  the  drive  motor  off  and  a  constant  torque  applied. 
The  same  behavior  persisted,  leading  to  the  belief  that  periodic 
voltage  surges  were  affecting  the  6940  multiprogrammer  unit. 
During  this  period  of  testing,  the  University  physical  plant  was 
shutting  down  at  night  all  nonessential  electrical  items,  such  as 
circulation  fans  in  the  ventilation  system,  in  the  building  in 
which  the  test  apparatus  was  located  to  conserve  energy.  Full 
power  was  brought  back  up  early  in  the  morning  and  it  was  thought 
that  this  may  have  been  the  cause  of  the  problem.  The  6940  unit 
was  removed  from  the  system  and  sent  to  be  serviced.  Nothing 
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could  be  found  defective  with  the  unit  so  it  was  again  installed 
and  testing  continued.  Further  tests  after  reassembly  showed  no 
abnormal  behavior.  No  satisfactory  answer  has  been  found  to 
explain  the  erratic  operation  discussed.  It  is  thought  that  the 
negative  torques  for  Test  64  are  caused  by  the  fact  that  the  out¬ 
put  is  subjected  to  some  error  and  the  voltage  surges  discussed. 

Test  65  was  run  with  a  carbon  seal  having  an  initial  conver¬ 
gent  taper  of  613  um/m.  The  torque  is  still  quite  low  at  0.68 
N*m  but  the  leakage  dropped  considerably  to  18  cm3/min.  This 
test  was  run  for  64  hours  and  displayed  the  jumps  in  torque  and 
seal  face  temperature  as  discussed  earlier. 

After  reinstalling  the  6940  multiprogrammer  unit.  Test  66 
was  run  for  118  hours.  The  initial  taper  was  48  um/m  convergent 
(nearly  initially  parallel  faces).  Test  66  in  the  Appendix  shows 
that  the  torque  increased  a  little  more  compared  to  the  previous 
test,  up  to  2.83  N-m.  The  leakage  has  again  dropped  consider¬ 
ably,  down  to  4  cm3/min.  The  operation  is  quite  smooth  except 
for  the  slightly  erratic  leakage  at  the  beginning. 

It  was  thought  that  a  parallel  face  operating  condition 
could  be  achieved  if  the  test  ran  long  enough.  However,  because 
of  thermal  rotation,  hydrostatic  effects  were  predominant  allow¬ 
ing  sufficient  lift  off  so  wear  was  quite  small.  After  118  hours 
torque  and  leakage  had  not  changed  enough  to  expect  parallel  face 
conditions  to  be  established  in  any  reasonable  length  of  time,  so 
the  test  was  ended. 

In  the  next  test  (Test  67,  Appendix),  the  carbon  seal  had  a 
divergent  taper  of  -320  um/m  lapped  into  it.  At  start  up,  the 
torque  was  measured  at  approximately  6.2  N*m  (54.8  in. -lb)  and 
then  decreased  rapidly  to  4.86  N*m  (43  in. -lb)  one  hour  later. 

The  torque  then  approached  a  constant  4.12  N*m  (36.5  in. -lb). 
Leakage  started  out  at  3  cm3/min  and  slowly  reduced  to  0.8 
cm3/min  after  32  hours  of  operation. 
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It  should  be  noted  here  that  the  delay  In  leakage  at  the 
start  of  each  test,  as  shown  In  the  Appendix  figures.  Is  caused 
by  the  fact  that  It  takes  about  20  cm3  of  leakage  to  fill  the 
inner  cavities  of  the  waviness  cylinder  and  leakage  catcher 
before  leakage  Is  able  to  flow  into  the  leakage  measuring  device. 

An  even  higher  divergent  taper  was  used  on  the  seal  for  Test 
68.  The  Initial  slope  was  -750  um/m.  Test  68  in  the  Appendix 
shows  that  the  initial  torque  starts  out  at  approximately  8.8  N»m 
(77.8  in. -lb)  and  drops  off  quite  fast  to  7.63  N*m  (67.5  in. -lb) 
after  one  hour.  Initial  leakage  is  about  2  cm3/min  and  slowly 
decreases  to  0.043  cm3/min  after  65  hours  of  operation.  The 
sharp  drop  in  torque  and  temperature  readings  at  approximately  19 
hours  was  due  to  a  momentary  (less  than  one  second)  shutdown  of 
the  test  apparatus.  The  motor  drive  was  not  allowed  to  stop 
before  power  was  reapplied.  It  appears  that  the  seal  was  able  to 
operate  in  two  stable  conditions  as  evidenced  by  the  torque  prior 
to  and  after  the  shutdown.  The  slopes  of  the  torque  and  tempera¬ 
ture  before  and  after  shutdown  are  quite  similar  although  leakage 
did  not  seem  to  be  affected.  The  final  torque  averaged  out  to 
7.71  N*m  (68.2  in. -lb)  after  65  hours. 

Test  75  was  run  with  a  divergent  taper  of  -1110  urn/ in.  The 
starting  trend  is  the  same  as  that  of  the  previous  tests.  The 
torque  started  out  at  approximately  11.8  N»m  (100  in. -lb)  and 
decreased  to  9.49  N*m  (84  in. -lb)  after  one  hour.  Operational 
trends  appeared  to  be  the  same  as  Tests  66  through  68.  However, 
after  four  hours  the  performance  was  quite  different  from  previ¬ 
ous  tests.  Apparently,  the  seal  had  reached  a  parallel  face  con¬ 
dition.  Leakage  had  dropped  off  very  fast  with  no  measurable 
leakage  after  eight  hours.  The  sudden  changes  in  torque  may  be 
caused  by  the  faces  drying  out,  the  result  of  the  fluid  being 
sealed  off  at  the  outside  diameter.  Since  dry  carbon  has  a 
higher  friction  coefficient  than  wet  carbon,  the  increased  fric¬ 
tion  results  in  a  sharp  increase  in  torque  and  temperature.  As 
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the  temperature  increases,  additional  thermal  rotation  occurs 
such  that  some  fluid  reenters  the  gap  and  friction  is  reduced 
[9].  The  final  corrected  torque  averaged  at  11.92  N»m  (105.5 
in. -lb)  after  71.5  hours  of  operation. 

Test  76  was  a  repeat  of  Test  68  so  as  to  obtain  additional 
wear  data.  The  initial  taper  was  -778  um/m.  Test  76  in  the 
Appendix  shows  the  initial  torque  at  6.89  N*m  (61  in. -lb)  and 
then  decreasing.  After  about  one  hour  the  torque  began  to  rise 
and  become  erratic.  Apparently,  a  parallel  face  condition  had 
developed  quite  quickly.  The  face  temperature  was  very  high  due 
to  the  fact  that  there  was  no  leakage  initially.  A  shutdown 
occurred  at  approximately  16  hours  into  the  test  due  to  excessive 
face  temperature.  The  maximum  level  was  reset  and  the  test  con¬ 
tinued. 

At  about  17  hours  leakage  began,  which  in  turn  reduced  the 
face  temperature  and  torque  dramatically.  The  amount  of  this 
initial  leakage  was  very  small  (approximately  0.1  cm3/min)  and 
shows  how  sensitive  torque  and  face  temperature  are  to  this  quan¬ 
tity.  As  the  leakagae  decreased,  the  torque  and  face  temperature 
once  again  began  to  increase.  The  performance  steadied  after  36 
hours.  The  final  corrected  torque  was  10.93  N*m  (96.74  in. -lb) 
and  leakage  was  0.00719  cm3/min.  These  values  were  averaged  over 
the  last  24  hours  of  operation. 

Effects  of  Pressure  Caused  Rotation 

After  Test  66  was  run,  a  separate  test  was  performed  to 
check  the  zero  moment  design,  i.e.,  a  check  to  insure  little  or 
no  rotation  of  the  seal  faces  occurred  due  to  sealed  pressure. 

Two  strain  gauges  were  mounted  on  the  seal  ring  as  shown  in  Fig¬ 
ure  2-9.  With  the  strain  gauges  in  these  locations,  a  divergent 
rotation  would  put  the  gauge  on  the  left  into  tension  and  the 
gauge  on  the  right  into  compression.  Calibration  was  accom¬ 
plished  by  the  use  of  the  pressure  type  radial  tapering  device. 
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Strain  Gauge 


Figure  2-9.  Carbon  Seal  with  Strain  Gages. 


The  device  was  inserted  into  the  seal  ring  toward  the  face.  A 
spacer  was  used  to  control  the  location  of  the  tapering  device 
relative  to  the  face.  After  the  tapering  device  was  pressurized, 
a  radial  surface  profile  was  made  across  both  faces  to  record  the 
face  rotation;  then  the  resultant  strain  was  measured  from  a 
strain  indicator.  The  results  of  calibration  are  shown  in  Figure 
2-10.  The  seal  was  then  mounted  on  the  seal  test  apparatus  and 
the  strain  gauge  leads  were  run  through  a  small  brass  tube  which 
exited  from  the  left-hand  end  plate  through  a  pressure  fitting. 
The  gauge  wires  were  epoxied  at  the  tube  opening  inside  the  pres¬ 
sure  vessel  to  prevent  leakage.  The  pressure  vessel  was  in¬ 
stalled  and  the  system  was  pressurized  with  water  to  3.45  MPa 
(500  psi).  The  results  of  the  experimentation  are  given  in  Table 
2-2.  The  average  tilt  was  103  um/m  divergent.  This  is  consid¬ 
ered  within  acceptable  limits  since  such  an  amount  would  affect 
the  performance  curves  very  little. 

Surface  Roughness 

All  mixed  friction  theoretical  models  to  date  are  based  on 
the  assumption  of  Gaussian  or  similar  surface  roughness.  It  is 
now  recognized  that  this  may  lead  to  considerable  error  in  the 
prediction  of  leakage  because  carbon  surfaces  are  not  Gaussian  by 
any  means. 

To  better  understand  the  nature  of  carbon  surfaces,  a  series 
of  measurements  was  made  as  a  part  of  the  radial  taper  test 
series.  Table  2-3  contains  the  surface  statistics  before  and 
after  the  test  in  both  the  radial  and  tangential  directions. 

These  data  were  obtained  using  a  surface  analyzer.  The  initial 
radial  and  tangential  a  are  quite  close,  leading  to  the  conclu¬ 
sion  that  the  inital  lapped  surface  characteristics  are  indepen¬ 
dent  of  direction.  The  final  values  show  a  considerable  reduc¬ 
tion  in  the  standard  deviation  in  both  the  radial  and  tangential 
directions. 
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Table  2-3.  Radial  Taper  Tests  Carbon  Surface  Statistics 


Test 

No. 

Initial  a 
Radial 

un  (yin.) 
Tangential 

Final  a 
Radial 

un  (yin.) 
Tangential 

64 

0.79  (31.3) 

0.83  (32.5) 

0.79  (31.3) 

0.70  (27.5) 

65 

0.70  (27.5) 

0.70  (27.5) 

0.64  (25.0) 

0.57  (22.5) 

66 

0.70  (27.5) 

0.79  (31.3) 

0.57  (22.5) 

0.51  (20.0) 

67 

0.73  (28.8) 

0.79  (31.3) 

0.70  (27.5) 

0.70  (27.5) 

68 

0.67  (26.3) 

0.76  (30.0) 

0.41  (16.3) 

0.38  (15.0) 

75 

0.73  (28.8) 

0.76  (30.0) 

0.35  (13.75) 

0.38  (15.0) 

Avg. 

0.72  (28.4)* 

0.77  (30.4) 

0.58  (22.7) 

0.54  (21.3) 
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With  these  values  and  using  the  mathematical  models  men¬ 
tioned,  one  would  expect  a  significant  leakage  in  the  parallel 
face  tests.  This  was  not  the  case  (see  Tests  75  and  76  and  all 
B  *  1  parallel  face  tests)  [8].  To  further  examine  this  ques¬ 
tion,  a  histogram  was  generated  for  the  carbon  surface  of  Test  68 
as  described  previously  [8]  and  is  shown  in  Figure  2-11.  An  ex¬ 
panded  surface  profile  is  shown  in  Figure  2-12.  The  histogram  is 
characteristic  of  a  worn  carbon  surface  having  a  series  of  val¬ 
leys  with  the  higher  asperities  being  chopped  off.  Figure  2-12 
also  shows  this  surface  characteristic. 

The  question  that  has  arisen  is  whether  or  not  the  valleys, 
as  seen  in  Figure  2-12,  are  connected  or  are  they  simply  isolated 
pockets.  If  they  are  connected,  the  measured  standard  deviation 
would  be  justified,  but  if  they  are  not  connected,  the  surface 
should  be  characterized  for  purposes  of  calculating  leakage  by 
the  much  smaller  standard  deviation  of  the  relatively  "flat" 
peaks.  For  the  sake  of  future  discussion,  the  idea  that  while 
the  valleys  contribute  to  the  measured  roughness,  they  do  not 
significantly  influence  flow  is  called  the  "pocket  theory." 

Proving  the  pocket  theory  by  use  of  the  trace  shown  in  Fig¬ 
ure  2-12  is  not  possible.  A  three-dimensional  view  is  needed. 

It  was  decided  that  an  electron  beam  microscope  scan  of  the  car¬ 
bon  surface  would  be  useful.  This  requires  the  destruction  of  a 
carbon  ring.  Seal  11,  which  was  used  for  Test  68,  was  chosen 
since  a  final  trace  of  the  surface  showed  that  the  inside  radius 
had  no  evidence  of  wear  while  the  remainder  of  the  surface  showed 
some  wear.  This  provided  the  opportunity  to  see  both  worn  and 
unworn  surface  characteristics  on  a  single  seal. 

Figure  2-13  is  the  image  of  the  unworn  section  viewed  radi¬ 
ally  from  inside  to  outside  at  an  angle  of  10*  from  the  horizon¬ 
tal.  The  dark  area  in  the  foreground  is  the  vertical  edge  at  the 
inside  radius  of  the  ring.  This  section  shows  the  inherently 
rough  nature  of  the  surface  due  to  lapping.  Figure  2-14,  on  the 
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Tangential  Roughness  Distribution  -  Test  68 
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Electron  Beam  Microscope  Image  of  Unworn  Section  -  Test  68 


Figure  2-14.  Electron  Beam  Microscope  Image  of  Worn  Section  -  Test  68. 


other  hand,  shows  the  worn  section  of  the  surface.  Again,  this 
is  viewed  radially  from  inside  to  outside  at  an  angle  of  10*  from 
the  horizontal.  Both  scans  are  magnified  1800  times.  Figure  2-14 
confirms  the  pocket  theory.  These  pockets  are  not  conducive  to 
fluid  flow  in  the  radial  direction  because  liquid  becomes  trapped 
in  them.  Therefore,  the  actual  surface  should  be  characterized 
from  a  leakage  standpoint  by  a  smaller  standard  deviation  than 
actually  measured. 

Radial  Surface  Profiles 

Figure  2-15  shows  the  initial  and  final  radial  traces  for 
Tests  64  and  65.  Looking  at  these  two  tests  and  the  correspond¬ 
ing  initial  and  final  surface  roughness  statistics  in  Table  2-3 
indicates  clearly  that  no  measurable  wear  had  taken  place. 
Although  the  final  slopes  are  not  exactly  the  same  as  the  ini¬ 
tial,  they  are  essentially  unchanged,  since  measurement  of  the 
slopes  is  no  better  than  ±50  um/m.  This  result  explains  why  the 
performance  curves  changed  very  little  during  testing. 

Figure  2-16  also  shows  no  measurable  wear  or  significant 
change  in  surface  roughness  for  both  Tests  66  and  67.  The  ther¬ 
mal  rotation  was  apparently  large  enough  in  Test  67  to  produce  a 
net  convergent  taper  during  operation  which  was  sufficient  to 
give  enough  hydrostatic  load  support  so  that  the  minimum  film 
thickness  was  greater  than  the  surface  roughness  and  hence  no 
wear.  The  initial  and  final  tapers  are  again  comparable. 

Tests  68  and  75,  Figure  2-17,  do  show  measurable  wear.  This 
is  also  seen  in  Table  2-3  where  the  surface  roughness  decreases 
from  the  beginning  of  the  test  to  the  end.  The  seal  in  Test  68 
has  not  completely  worn  to  a  flat  face  condition.  The  inside 
radius  shows  little  wear  on  this  trace.  The  asperities  on  the 
inside  were  still  quite  distinct  while  the  remainder  of  the  face 
shows  that  these  higher  asperities  had  been  chopped  off.  The 
electron  microscope  images.  Figures  2-13  and  2-14,  were  made  from 


this  surface.  The  final  trace  in  Test  68  does  not  represent  the 
equilibrium  thermal  rotation,  but  Is  probably  close.  As  the  seal 
continues  to  wear  in,  the  friction  will  increase  and  so  will  the 
thermal  rotation.  Therefore,  it  would  be  expected  that  the  final 
tilt  would  be  slightly  more  divergent.  The  flat  portions  on  the 
outside  edge  were  probably  caused  by  initial  start-up  and  finding 
the  zero  torque  level. 

The  final  trace  of  Test  75  shows  a  concave  shape.  Determin¬ 
ing  the  taper  is  quite  difficult  in  this  case,  so  the  profile  of 
the  mating  tungsten  carbide  ring  is  shown  in  Figure  2-18.  The 
tungsten  carbide  ring  has  deep  circumferential  grooves  that  cor¬ 
respond  to  the  exaggerated  peaks  at  the  inner  and  outer  radii  of 
the  carbon  ring.  Figure  2-18  shows  how  these  sections  match  up. 
By  taking  a  reference  point  at  the  outside  radius  as  being  zero 
and  finding  the  difference  between  heights  of  the  two  surfaces 
from  that  reference,  a  resultant  profile  is  obtained  and  is  pre¬ 
sented  in  Figure  2-19.  This  trace  still  poses  some  difficulty  in 
determining  a  taper  so  a  least  squares  fit  was  made.  The  final 
average  taper  was  -710  um/m  and  is  at  best  an  approximation  under 
the  circumstances.  The  reason  for  this  concave  shape  Is  due  most 
likely  to  the  fact  that  the  seal  face  temperature  was  quite  high 
locally  and  local  thermal  deformations  were  superimposed  over  the 
thermally  caused  tilt.  Tests  68  and  75  together  show  that  the 
thermal  rotation  at  the  parallel  face  operating  condition  is  from 
700-800  um/m.  This  will  be  compared  to  theory  later. 

Test  76,  Figure  2-20,  is  a  repeat  of  Test  68  to  get  addi¬ 
tional  wear  data.  Once  it  was  determined  from  Tests  68  and  75 
that  the  equilibrium  thermal  rotation  for  parallel  faces  is  about 
-750  um/m.  Test  76  was  run  at  an  initial  taper  of  -778  um/m. 
Figure  2-20  shows  a  thermal  taper  which  is  higher  than  the  re¬ 
sults  above.  For  some  unexplained  reason.  Test  76  ended  up  oper¬ 
ating  at  a  higher  torque  than  Tests  68  or  75,  thus  resulting  in 
the  higher  thermal  taper. 


19.  Resultant  Profile  -  Test  75. 


Wear  measurements  were  made  for  both  Tests  75  and  76.  The 
wear  values  are  7.04  and  10.67  un  per  100  hours,  respectively. 
Both  of  these  values  are  lower  than  obtained  for  the  B  =  1.0 
tests  as  shown  in  the  previous  annual  report  [8]. 

High  Temperature  Tests 

A  series  of  twenty-nine  tests  were  conducted  in  water  with 
nominally  flat  face  seals  at  elevated  temperatures  up  to  190*C 
for  both  B  *  0.75  and  B  =  1.00.  The  purpose  of  these  tests  was 
to  provide  data  to  compare  to  results  predicted  by  a  two-phase 
model  [11]  which  was  developed  in  part  during  this  research  pro¬ 
gram.  In  essence  the  model  predicts  that  seal  friction  should 
decrease  with  increasing  temperatures  because  more  of  the  load  is 
carried  by  fluid  pressure  as  the  phase  change  boundary  moves 
across  the  seal  face.  This  trend  will  continue  up  to  some  criti¬ 
cal  temperature  after  which  the  friction  will  again  increase. 
Theoretically,  leakage  will  suddenly  increase  at  this  critical 
point  and  seal  puffing  may  occur. 

Several  modifications  of  the  test  apparatus  were  required  to 
make  these  tests.  The  first  of  these  was  to  provide  a  means  to 
maintain  the  high  temperatures  required.  After  investigating  the 
heat  loss  through  the  pressure  vessel  and  the  space  available 
inside  the  pressure  vessel,  it  was  determined  that  a  tubular 
heating  element  with  a  capability  of  2  kW  was  adequate.  It  was 
bent  and  placed  inside  the  pressure  vessel  to  heat  the  system 
(Figure  2-21).  The  water  temperature  is  sensed  by  a  thermocouple 
and  monitored  by  computer.  Whenever  the  water  temperature  goes 
above  the  set  temperature,  the  control  program  opens  a  relay 
which  stops  the  electrical  heating. 

It  was  discovered  early  in  this  test  series  that  both  heat¬ 
ing  and  cooling  were  necessary  to  be  able  to  control  the  tempera¬ 
ture  under  all  operating  conditions.  At  the  lower  test  tempera¬ 
tures  the  cooling  system  described  previously  [8]  was  used. 
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The  cooling  control  temperature  was  set  a  few  degrees  above  the 
heating  temperature.  Since  the  water  temperature  was  higher  than 
the  temperature  limit  of  the  circulation  pump  for  some  of  the 
high  temperature  tests,  the  circulation  pump  was  by-passed  to 
avoid  damage.  A  small  re-circulation  vane  was  attached  on  the 
left-hand  end  plate  of  the  pressure  vessel  to  increase  circula¬ 
tion  and  cooling.  Experimental  results  show  that  temperature 
control  is  in  the  range  of  ±3*C  for  most  of  the  tests. 

When  conducting  a  test  at  a  high  temperature,  vapor  leakage 
would  be  expected  when  two-phase  operation  is  occurring.  After 
the  vapor  passes  the  seal  face,  it  will  condense  on  cooler  seal 
parts,  and  the  observed  leakage  will  be  in  the  form  of  a  liquid. 

A  newly  designed  leakage  catcher  (Figure  2-22)  with  an  0-ring  at 
the  right-hand  end  of  the  leakage  sleeve  was  used  during  the 
test.  The  vapor  leakage  is  sealed  and  condensed  to  a  liquid 
within  the  sleeve.  Then,  the  entire  leakage  is  collected  and 
measured  by  the  leakage  measuring  device. 

Based  on  the  experimental  evidence,  the  minimum  measurable 
leakage  for  the  measuring  device  is  about  0.05  cm3/min.  Below 
such  rates  it  is  expected  that  water  evaporates  from  the  leakage 
counter  before  it  can  be  registered.  Since  no  leakage  was  de¬ 
tected  for  all  the  tests  after  the  seals  were  worn  in,  leakage 
rates  below  this  value  must  have  occurred. 

A  Type  J  thermocouple  was  used  to  measure  the  seal  face 
temperature.  The  thermocouple  is  placed  at  a  specific  position 
beneath  the  seal  face  (Figure  2-23).  The  thermocouple  is  read 
several  times  per  minute,  converted  to  temperature,  averaged, 
plotted,  and  the  data  stored  on  tape.  There  are  three  sources  of 
error  for  temperature  measurement:  1)  the  accuracy  of  the  thermo¬ 
couple  (±2.2*C);  2)  the  accuracy  of  the  low-level  A/D  card 
(±0.7*C);  and  3)  the  error  caused  by  the  thermocouple  position 
error.  Based  upon  the  theoretical  temperature  distribution, 
there  is  a  temperature  gradient  of  about  30* C/cm  near  the  end  of 
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thermocouple.  The  thermocouple  position  error  is  about  ±0.025 
cm.  This  will  cause  an  error  of  ±0.75*C.  Therefore,  the  total 
absolute  error  for  seal  face  temperature  measurement  could  be  as 
much  as  ±3.7*C.  The  comparative  error  from  one  test  to  another 
and  within  one  test  is  much  smaller. 

In  all  of  the  experiments  performed,  the  outside  pressure 
and  the  drive  speed  were  adjusted  to  1.72  MPa  absolute  (250  psia) 
and  1800  rpm,  respectively.  Ordinarily,  the  seal  was  started 
with  the  desired  pressure  and  accelerated  immediately  to  the 
operating  speed.  The  system  was  heated  gradually  up  to  the 
desired  operating  temperature  within  about  one  hour.  Test  dura¬ 
tion  varied  from  1  hour  to  89  hours.  The  seal  was  lapped  only 
once  before  installation  and  run  through  several  tests.  The  ini¬ 
tial  surface  parameters  of  test  seals  are  shown  in  Table  2-4. 

Immediately  after  starting  a  new  or  relapped  seal,  there  was 
a  certain  period  characterized  by  relatively  low  torque  and  seal 
face  temperature  (see  Appendix)  during  which  leakage  was  de¬ 
tected.  This  behavior  agrees  with  the  theory  presented  in  refer¬ 
ence  [9]  which  suggests  that  during  initial  operation,  thermal 
coning  produces  a  convergent  film  shape  and  reduces  friction.  As 
the  face  wears  flat,  the  friction  increases.  Within  10  hours  or 
less,  conditions  stabilized  and  the  seals  appeared  to  have 
run-in. 

No  measurable  leakage  was  detected  after  the  run-in  period. 
After  the  seals  had  worn  in,  a  fluctuating  torque  and  seal  face 
temperature  were  recorded.  The  mean  value  of  torque  and  seal 
face  temperature  stabilized  at  a  certain  level.  At  this  point, 
it  was  assumed  that  the  seal  faces  had  been  worn  parallel. 

Tests  were  conducted  on  four  seals  with  operating  tempera¬ 
ture  varying  from  37.8*C  (100*F)  to  190. 6*C  (375’F)  as  shown  in 
Tables  2-5  and  2-6.  Most  of  the  data  are  averages  over  the  last 
8  hours  of  the  test.  The  Appendix  contains  all  the  detailed  per¬ 
formance  data  for  these  tests. 
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Table  2-4.  Initial  Surface  Data,  High  Temperature  Tests 


Seal  No. 

Radial  Taper 

Tangential 
CLA  (urn) 

Radial 

CLA  (wn) 

5 

Flat 

0.51 

0.58 

7 

Flat 

0.43 

0.66 

9 

Flat 

0.43 

0.28 

10 

-570.10-6 

0.58 

0.61 

10* 

Flat 

0.46 

0.41 

♦Lapped  after  Test  73. 
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Table  2-6.  High  Temperature  Test  Experimental  Results, 
B  »  0.75,  Ph2q  *  1-72  MPa,  1800  rpm 


Test 

No. 

Seal 

No. 

Test 

Duration 

(h) 

H 

Average 

Torque 

(N-m) 

Average 
Face  Temp 
CC) 

Average 

Leakage 

(cm3/min) 

69 

10 

24.1 

37.8 

4.6 

44.1 

mm 

70 

10 

24.1 

93.3 

6.9 

98.1 

BB 

71 

10 

24.0 

148.9 

6.7 

153.2 

■H; 

72 

10 

89.2 

37.8 

1.7 

41.1 

0.0 

73 

10 

5.0 

37.8 

to 

190.6 

78 

10* 

37.1 

37.8 

7.6 

47.1 

0.0 

79 

10 

12.2 

65.6 

5.3 

71.0 

0.0 

80 

10 

10.0 

93.3 

7.4 

100.8 

0.0 

81 

10 

12.0 

121.1 

6.5 

126.7 

0.0 

82 

10 

11.7 

176.7 

4.2 

180.4 

0.0 

83 

10 

3.8 

148.9 

4.0 

152.8 

0.1 

84 

10 

2.1 

190.6 

10.1 

199.7 

0.0 

♦Lapped  flat  after  Test  73. 
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Tests  58  and  59  are  short  because  of  a  shut  down  due  to 
overload  of  the  drive  motor  immediately  after  the  operating  tem¬ 
perature  reached  176. 7*C  (350*F).  Three  short  period  tests 
(Tests  62,  63,  and  73)  were  run  for  the  purpose  of  observing  the 
effects  of  operating  temperature  on  seal  performance  on  a  contin¬ 
uous  basis.  Tests  62  and  63  were  shut  down  by  control  program 
due  to  overload  of  the  drive  motor  at  T*,  *  176. 7*C  (350*F). 

Test  73  was  shut  down  manually  when  puffing  occurred  at  T„  = 

191 *C . 

In  order  to  verify  the  parallel  face  operation,  an  examina¬ 
tion  of  the  torque  curves  was  made  after  all  tests  had  been  com¬ 
pleted.  Results  show  that  most  of  the  torque  curves,  except  for 
Tests  69,  72,  and  79  are  characterized  by  a  wide  variation  in 
torque  level  (see  Appendix),  which  is  characteristic  of  a  paral¬ 
lel  face  operation.  The  torque  curves  of  Tests  59,  72,  and  79 
were  quite  steady  and  the  torque  levels  were  much  lower  compared 
to  the  others.  It  is  thought  that  these  seals  were  operating 
with  a  convergent  taper,  and  this  resulted  in  a  sufficient  hydro¬ 
static  load  support  which  prevented  wear-in.  Thus,  they  are 
excluded  from  the  data  which  will  be  used  for  later  comparison  to 
theory. 

Figures  2-24  and  2-25  show  the  results  of  the  continuous 
Tests  62  and  73.  In  these  tests  temperatures  were  gradually 
increased  with  time.  The  data  show  the  decreasing  and  then 
increasing  trend  in  aT  (or  equivalently  torque)  as  discussed 
earlier. 

Surface  profiles  were  made  before  and  after  each  test  series 
using  the  surface  analyzer.  The  divergent  radial  taper  caused  by 
the  thermal  coning  effect,  shown  in  Figure  2-26,  agrees  with  the 
theory  presented  in  reference  [9].  But,  the  convergent  profile 
shown  typically  by  Figure  2-27  is  in  disagreement  with  theory. 
Close  examination  of  the  operational  history  of  each  test  seal 
shows  that  the  final  operating  temperature  for  the  convergent 
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Figure  2-2 5.  High  Temperature  Continuous  Test,  B  =  0.75. 
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Final 


Figure  2-26.  Radial  Surface  Profile  -  Seal  #5 
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Figure  2-27.  Radial  Surface  Profile  -  Seal  #7, 


seal  was  above  176. 7*C  (350*F),  whereas  for  the  divergent  seal  it 
was  100*F  (30*C).  It  was  thought  that  at  high  operating  tempera¬ 
tures,  the  thermal  expansion  of  the  waviness  ring  might  cause 
interference  with  the  seal  and  cause  it  to  become  divergent  dur¬ 
ing  operating.  To  investigate  this,  the  inside  dianeter  of  a 
carbon  ring  (Seal  10*)  was  increased  by  0.5  mm  and  run  through 
several  tests.  The  final  trace  still  indicates  a  convergent  pro¬ 
file.  This  means  that  there  must  be  another  kind  of  thermal  dis¬ 
tortion  taking  place  for  high  temperature  operation.  Further  in¬ 
vestigation  is  needed  to  understand  how  the  carbon  ring  deforms 
in  the  high  temperature  environment. 

Although  this  behavior  may  have  influenced  the  test  results, 
in  most  tests  the  period  of  operation  was  long  enough  that  the 
faces  became  parallel  in  spite  of  the  additional  apparent 
thermally-caused  rotation. 


Long-Term  Tests 

Wear  tests  results  on  the  wavy  seals  reported  in  the  previ¬ 
ous  report  [8]  show  some  inconsistency  and  in  some  cases  are 
clearly  beyond  the  capability  of  measurement  techniques.  To  get 
better  wear  data,  it  was  decided  to  run  a  series  of  long-term 
wear  tests  of  500  hours  each.  Table  2-7  shows  the  results  of  two 
such  tests.  Test  45  was  run  using  the  original  waviness  drive 
concept  with  sliding  0-rings,  whereas  Test  110  was  run  using  the 
new  waviness  generator  described  previously.  Comparing  the  ac¬ 
tual  test  data  between  Test  45  [8]  and  Test  110  (see  Appendix) 
shows  that  much  more  uniform  performance  was  obtained  using  the 
new  wav i ness  generator  than  the  previous  method.  The  leakage 
oscillations  in  Test  110  are  much  smaller.  Recalling  that  the 
oscillations  are  caused  by  the  superpositions  of  the  second  har¬ 
monic  drive  force  caused  waviness  over  the  three  intentional 
waves,  the  oscillation  in  Test  110  is  smaller  because  torque  on 
the  drive  lugs  is  small  since  there  is  no  friction  torque  on  the 
seal  ring  caused  by  sliding  gas  pocket  0-ring  friction. 


Table  2-7.  Long-Term  Wear  Test  Results 
1800  rpm,  pQ  =  3.45  MPa,  pg  a  6.90  MPa,  8  *  1,  n  »  3 


Test  No. 

Duration 

00 

Torque* 

N*m 

Leakage* 

(ml/min) 

Face  Temp 
(°C) 

Avg.  Carbon 
Wear 
( wn) 

45 

507 

3.39 

6.8 

32.2 

0.31 

110 

500 

1.28 

1.54 

39.3 

-0.84 

♦Average  o^er  last  9.47  hours  of  test. 


Average  leakage  for  the  final  9.47  hours  of  Test  110  is 
lower  than  for  Test  45.  It  is  thought  that  leakage  is  smaller 
because  there  are  no  large  oscillations  as  in  Test  45.  Test  45 
in  fact  shows  that  the  oscillation  amplitude  increases  toward  the 
end  of  the  test  as  0-ring  friction  presumably  becomes  higher. 

The  opposite  is  true  of  Test  110  (see  Appendix).  After  the  first 
24  hours  the  amplitude  of  the  leakage  oscillations  decrease  as 
the  seal  wears  in.  Likewise,  torque  in  Test  110  is  lower  than  in 
Test  45. 

Table  2-7  shows  the  wear  rate  for  Test  45  to  be  very  small, 
0.3  pm/500  hours.  Test  110  shows  a  negative  wear  rate  meaning 
that  wear  is  less  than  the  distortion  of  the  nose  of  the  seal. 
Comparing  the  lowest  wear  rate  from  the  parallel  face  tests  of 
the  previous  report  [8]  to  the  worst  case  above  (Test  45)  shows 
that  the  wear  rate  for  the  wavy  seal  is  a  factor  of  at  least  200 
times  lower  than  for  the  parallel  face  seal.  At  such  a  low  wear 
rate  a  wavy  seal  could  be  designed  such  that  a  ten  year  face  life 
in  continuous  operation  could  be  anticipated. 

Offset  Land  Tests 

As  will  be  discussed  later,  it  can  be  shown  theoretically 
that  a  flat  face  seal  constructed  with  an  offset  will  have  a 
lower  wear  rate  than  a  parallel  face  seal.  This  seal  can  be  made 
to  operate  theoretically  with  a  low  leakage  rate  for  an  indefin¬ 
ite  period  of  time  by  in  essence  machining  the  offset  into  the 
soft  face  on  a  continuous  basis.  This  concept  is  shown  in  Figure 
2-28.  The  hard  face  operates  hydrostatically  as  if  it  were  a 
flat  face.  The  only  purpose  of  the  protrusion  is  to  machine  the 
offset  into  the  mating  soft  face  as  shown. 

To  produce  such  a  seal  it  is  necessary  to  provide  the  pro¬ 
trusion  as  shown  while  the  remainder  of  the  face  is  left  undis¬ 
turbed.  The  height  of  the  protrusion  is  of  the  order  of  1-5  pm 
depending  on  performance  desired.  As  a  first  attempt  an  electron 
beam  deposition  device  was  used  to  produce  the  protrusion  on  the 
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hard  face.  Tungsten  carbide  powder  was  vaporized  and  deposited 
on  the  face  using  a  mask  to  obtain  the  desired  pattern  of  protru- 
sions  similar  to  that  shown.  It  was  found  that  the  deposits 
would  not  bond  tightly  to  the  seal  face.  After  several  trials 
using  different  techniques,  this  method  was  abandoned. 

The  technique  which  was  found  to  work  was  to  use  the  diamond 
indentor  from  a  Rockwell  tester  under  load.  This  technique  pro¬ 
duces  a  protrusion  of  the  form  shown  by  the  surface  analyzer 
trace  of  Figure  2-29.  The  height  of  the  protrusion  can  be  con¬ 
trolled  by  lapping  with  diamond  compound  and  by  varying  the 
indentor  load.  The  protrusion  produced  in  this  way  is  actually  a 
round  crater  with  the  actual  protrusion  being  formed  by  the  walls 
of  the  crater.  This  type  of  protrusion  has  the  disadvantage  that 
in  nonductile  materials  the  crater  walls  may  have  a  tendency  to 
break  off  as  discussed  later. 

Using  this  technique,  seal  faces  were  prepared  with  A/B  = 

0.2  as  shown  in  Figure  2-28,  and  four  tests  were  run  as  shown  in 
Table  2-8.  The  Appendix  shows  the  actual  performance  data  for 
these  tests.  Test  105  was  started  up  with  no  special  procedure. 
The  running  torque  is,  in  fact,  higher  than  observed  on  previous 
tests  for  parallel  face  seals  and  the  theoretical  mechanism  was 
clearly  not  operating.  Test  106  was  started  at  400  rpm  under 
full  water  pressure  and  run  for  approximately  5  minutes.  Then 
speed  was  increased  to  1800  rpm.  Again  here,  torque  was  higher 
than  in  parallel  face  tests.  In  both  Tests  105  and  106  it 
appeared  that  the  protrusions  were  acting  to  grind  away  the  mat¬ 
ing  face  and,  in  fact,  causing  a  higher  torque  than  for  the  par-  ■ 
a  1 1 e 1  face  test.  Judging  by  the  low  leakage,  the  fluid  never  had 
the  opportunity  to  become  distributed  properly  across  the  face  to 
relieve  the  mechanical  load. 

Test  107  was  set  up  with  much  larger  protrusions.  It  was 
started  up  under  about  60  psi  and  run  at  400  rpm  for  nearly  one 
hour  after  which  speed  and  pressure  were  increased  to  test 


Figure  2-29.  Surface  Analyzer  Trace  of  Diamond 
Indentor-Produced  Protrusions. 
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values.  In  contrast  this  seal  clear ly  lifted  off.  Leakage  was 
quite  high  and  torque  was  extremely  low.  Test  108  was  run  in  an 
identical  manner  but  using  a  smaller  protrusion.  Leakage  became 
more  moderate  and  torque  was  still  quite  low.  Figure  2-30  shows 
the  final  offset  and  protrusion  for  Test  108,  and  it  can  be 
clearly  seen  how  the  offset  is  machined  into  the  carbon.  The 
data  for  Test  108  also  shows  that  wear  is  quite  low  as  predicted 
even  for  the  162  hour  test.  Table  2-8  also  shows  that  the  pro¬ 
trusions  wear  down. 

In  summary,  the  concept  appears  to  work.  However,  there  are 
certain  limitations  and  unknowns.  First,  from  Test  108,  the  pro¬ 
trusions  appear  to  wear  down.  It  is  not  known  whether  they  will 
wear  down  to  a  still  useful  level  or  simply  wear  away  over  a  long 
period  of  time.  Second,  observing  that  Tests  105,  106,  and  107 
have  the  same  final  W-C  offset,  there  may  be  serious  problems  in 
repeatedly  starting  up  such  a  seal  in  the  low-friction  mode.  The 
seal  can  clearly  start  also  in  a  high-friction  mode,  at  least  at 
the  lower  offset  values.  This  is  a  real  disadvantage  because 
only  at  low  offset  values  will  the  leakage  assume  a  low  value— 
which  is  one  of  the  theoretical  advantages  of  this  seal  over 
other  types  of  hydrostatic  seals  which  characteristically  have 
high  leakage.  Further  work  is  needed  to  see  what  is  required  to 
make  the  seal  start  up  in  the  proper  mode.  One  possibility  to  be 
considered  is  to  use  an  offset  taper  arrangement  which  has  a  sim¬ 
ilar  theoretical  performance  but  may  have  a  stronger  tendency  to 
start  properly. 

Speed  Effects 

Several  additional  tests  were  run  at  varying  speeds  for  both 
the  wavy  seal  and  the  flat  face  seal.  The  results  of  these  tests 
plus  similar  tests  run  during  the  previous  year  are  shown  in 
Tables  2-9  and  2-10.  In  Table  2-9  for  the  wavy  tests,  the  most 
striking  result  is  the  very  high  friction  coefficient  at  1.4  rpm 
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Table  2-9.  Wavy  Seal  Speed  Effect  at 
pQ  =■  3.45  MPa,  pg  =  6.90  MPa,  n  -  3,  B  -  1.0 


Test 

No. 

Speed 

RPM 

Test 

Duration 

00 

Avg.  Face 
Temp 

C  c) 

Avg. 

Torque 

N-M  (in. -lb) 

Avg. 

Leakage 

cnr/min. 

91 

1.4 

20.6 

37 

29.0 

(255) 

1.2 

103 

1.4 

48 

38 

30.7 

(272) 

1.3 

90 

400 

72.4 

38 

0.6 

(6) 

2.3 

102 

400 

47.8 

38 

0.7 

(6) 

1.9 

26 

500 

4 

40 

1.6 

(14) 

0.3 

33 

900 

20 

40 

1.8 

(16) 

0.2 

104 

900 

28.7 

38 

0.3 

(3) 

1.4 

89 

900 

50 

39 

2.2 

(19) 

0.6 

34 

1800 

100 

39 

1.4 

(12) 

1.4 

39 

1800 

100 

39 

1.6 

(14) 

3.2 

21 

1800 

100 

40 

1.5 

(13) 

1.1 

30 

2700 

46 

41 

3.2 

(28) 

4.5 

31 

3600 

6 

43 

5.0 

(44) 

1.8 
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Table  2-10.  Flat  Face  Speed  Effects, 
pQ  =  3.45  MPa,  B  =  1.0 


Test 

No. 

Speed 

PPM 

Test 

Duration 

(h) 

Avg.  Face 
Temp 

(*c) 

Avg. 

Torque 

N*M  (in. -lb) 

Avg. 

Leak  age 
cm^/min. 

92 

1.4 

13.1 

37 

18.6 

(164) 

0.0 

93 

400 

23.6 

41 

9.1 

(80) 

0.0 

94 

900 

25.1 

43 

8.6 

(76) 

0.0 

35 

1800 

100 

44 

8.4 

(74) 

0.0 

40 

1800 

i 

111 

43 

9.2 

(81) 

0.0 

95 

2700 

22 

47 

6.2 

(54) 

0.0 

96 

3600 

0.3 

50 

6.2 

(55) 

0.0 
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and  the  dramatic  reduction  in  torque  at  a  relatively  slow  speed 
of  400  rpm.  After  400  rpm,  friction  torque  generally  increases 
with  increasing  speed. 

Looking  at  the  flat  face  speed  effects  in  Table  2-10,  again 
torque  is  very  high  at  1.4  rpm,  drops  off  rapidly  at  400  rpm,  and 
then  continues  to  decrease  with  increasing  speed.  The  flat  face 
results  can  be  explained  in  two  ways.  Assuming  that  the  seal  has 
worn  to  a  parallel  face  configuration  in  all  cases,  then  there 
must  exist  a  reduction  in  mechanical  coefficient  of  friction  with 
increasing  speed  or  microasperity  hydrodynamic  lubrication  re¬ 
duces  the  effective  mechanical  coefficient  of  friction,  or  both. 
Unfortunately,  there  is  no  way  to  distinguish  between  the  two 
effects.  Alternatively,  part  of  the  decrease  might  be  attribut¬ 
able  to  an  increase  in  thermal  coning  with  increasing  speed  if 
the  faces  are  not  worn  parallel  in  each  instance.  However,  this 
effect  could  not  account  for  the  large  difference  between  1.4  and 
1800  rpm  because  it  is  known  that  the  1800  rpm  cases  have  been 
well  worn  in.  Thus,  this  large  difference  must  be  explained  by 
one  of  the  first  two  mechanisms.  A  better  distinction  between 
the  thermal  coning  mechanism  and  changing  friction  coefficient 
will  become  available  using  the  friction  and  wear  tester  dis¬ 
cussed  previously  because  thermal  coning  and  related  hydrostatic 
effects  can  be  eliminated. 

Now  considering  the  trends  in  the  waviness  tests  in  Table 
2-9,  the  sane  effect  as  in  the  flat  face  test  must  be  operating. 
However,  the  decrease  is  greater,  and  this  is  attributable  to  the 
greater  global  hydrodyanamic  effect  caused  by  the  waves.  The 
reason  for  the  subsequent  increase  can  be  explained  in  part  by 
the  fact  that  the  wave  changes  in  amplitude  with  speed  as  dis¬ 
cussed  later  in  Chapter  4.  It  is  expected  that  when  the  nine- 
wave,  redesigned  seal  is  tested,  this  confusing  behavior  will  be 
eliminated  and  the  wavy  seal  speed-dependent  behavior  as  caused 
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by  hydrodynamic  action  will  be  established  conclusively. 

Nine-Wave  Tests 

The  waviness  cylinder  under  the  seal  ring  was  modified  to 
have  18  pressure  pockets.  If  each  pocket  is  pressurized  equally, 
nine  waves  result  and  these  can  be  moved  in  the  original  fashion 
by  rotating  the  waviness  cylinder..  Thus,  Table  2-11  shows  four 
tests  conducted  using  nine  waves.  In  all  cases  friction  is  quite 
high  and  leakage  is  zero  or  quite  small.  This  indicates  that 
neither  significant  enhanced  hydrostatic  support  or  hydrodynamic 
load  support  occurred.  As  will  be  discussed  later,  analysis  pre¬ 
dicts  this  behavior  in  general.  It  occurs  because  the  resultant 
deflection  of  the  original  seal  ring  deflected  to  nine  waves 
takes  on  a  disadvantageous  shape  and  does  not,  in  fact,  mean  that 
a  nine-wave  shape  of  the  proper  deflection  characteristic  will 
behave  poorly.  This  case  serves  as  an  example  that  simply 
deflecting  a  seal  into  a  wavy  shape  will  not  necessarily  be  of 
advantage.  The  exact  nature  of  the  shape  that  creates  an  advan¬ 
tage  must  be  understood  and  properly  imposed. 
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Table  2-11.  Wavy  Seal  Results,  n  =  9, 
p0  =  3.45  MPa,  pg  =  6.90  MPa,  B  =  1.0 


Test 

No. 

Speed 

RPM 

Test 

Duration 

00 

Avg.  Face 
Temp 

r  c) 

Avg. 

Torque 

N*m  (in. -lb) 

Avg. 

Leakage 

cm-fymin. 

85 

1800 

24.5 

53 

13.2 

(117) 

0 

86 

900 

48.5 

40 

5.1 

(45) 

0 

87 

1800 

44.3 

46 

7.2 

(64) 

0 

88 

2700 

44.8 

53 

8.7 

(77) 

0.03 
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Chapter  3 

Seal  Theory  Developments 

Waviness  Drive 

It  became  clear  soon  after  the  moving  wavy  seal  was  first 
tested  that  moving  the  wave  by  moving  pressure  pockets  with  slid¬ 
ing  0-rings  represented  a  serious  limitation  on  what  Is  otherwise 
an  extremely  long-life  seal  design.  The  500-hour  test  (Test  45) 
shown  In  the  previous  report  [8]  Is  a  case  In  point.  Toward  the 
end  of  the  test,  leakage  amplitude  became  Urge  and  erratic. 

After  disassembly  the  0-rings  showed  considerable  wear,  and  it  is 
thought  that  0-ring  friction  became  high  toward  the  end  of  the 
test.  This  caused  the  second  harmonic  distortion  to  become  Urge 
giving  rise  to  a  larger  leakage  amplitude. 

Although  other  mechanical  means  of  moving  the  wave,  such  as 
pressure  loaded  rollers,  were  considered,  it  was  decided  that  the 
most  reliable  operation  could  be  obtained  if  there  were  no  Inter¬ 
nal  moving  parts  except  the  seal  faces  themselves.  There  are 
various  concepts  which  can  be  used  to  move  the  distortion  pattern 
with  no  internal  moving  parts.  The  method  selected  for  develop¬ 
ment  Is  similar  conceptually  to  that  of  a  three-phase  electric 
motor  where  three  voltages,  each  at  a  different  phase  angle,  are 
applied  to  three  sets  of  windings  giving  rise  to  a  rotating  mag¬ 
netic  field.  In  the  case  of  the  seal,  three  pressures,  each  dif¬ 
ferent  in  phase,  are  applied  to  three  sets  of  pressure  pockets 
causing  the  resulting  distortion  to  move. 

With  reference  to  Figure  3-1,  three  sets  of  pressure  pockets 
similar  to  those  described  previously  [8]  are  shown.  Each  set  of 
six  pockets  A,  8,  or  C  will  produce  three  waves  around  the  seal 
because  of  the  alternating  moment  produced  on  the  ring.  It  can 
be  shovwi  that  the  first  and  most  significant  term  of  the  Fourier 
series  approximation  to  the  moment  produced  by  each  set  of 
pockets  is  given  by  [8] 
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Figure  3-1.  Pressure  Pocket  Arrangement  for  Moving  Waviness  Generator. 


meA  3  PA  sin('2?7)  cos  n0  *  (3_1) 

ffl0B  3  PB  “rj"  Sin(^)  C0S(n9  ' t)  *  {3‘2) 

mec 3  pc  “F”  sin(iF:)  cos(ne  ‘  t)  *  (3_3) 

C  I 

where 

n  *  3  (3-4) 


for  the  example  in  Figure  3-1.  The  predominant  term  for  deflec 
tion  due  to  an  nth  harmonic  distribution  such  as  mQ  is 


rj(l  +  A) 

v  = - - - *■  m  . 

EJx(n2  -  l)2  9 

Therefore,  the  total  deflection  v  at  any  9  is 
r|(l  +  A) 

V  3  EJ  (n2  -'17  (m0A  +  m0B  +  m0c)  ’ 


(3-5) 


(3-6) 


r;(l  +  A)  dr . e  /nH  x 

£J  ‘n8  7^7  —  s1n(^7)  [»A  cos  "9  +  pB  cos  ("9  -  T ) 


/  4  IT 

+  p^  cosing  -  y- 


)]• 


(3-7) 


Now  assume  the  pressure  pA,  pg,  and  pc  are  functions  of  time 
as  follows: 


PA  *  pg  cos  nt  +  p0  , 

PB  =  pg  “*(*  +TL)  +  po  ’ 


(3-8) 

(3-9) 


i 
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(3-10) 


PC  *  pg  COs((Jt  +y!)+  P„  • 

Substitution  of  equations  (3-8),  (3-9),  and  (3-10)  into  (3-7) 
gives 

v  *  K  cos(ne  -  nt)  ,  (3-11) 

where 


A  similar  derivation  applies  for  the  face  tilt  <j>.  Thus,  equation 
(3-11)  represents  a  deflection  of  n  waves  which  moves  at  angular 
speed  fl  around  the  face  of  the  seal.  This  is  accomplished  by 
supplying  three  pressures,  varying  cosinusoidally  with  time  and 
spaced  120*  apart,  to  three  separate  sets  of  pockets.  The  theory 
equally  applies  to  any  number  of  waves  n,  only  the  number  of  poc¬ 
kets  need  be  changed.  Three  pressures  only  are  required. 

Pressures  of  the  form  of  equation  (3-8),  (3-9),  and  (3-10) 
can  be  obtained  by  several  means  including  the  waviness  genera¬ 
tor  described  in  detail  in  Chapter  2.  It  is  to  be  noted  that 
each  pressure  must  have  the  same  amplitude  pg  and  the  same  mean 
value  p0.  If  these  are  not  equal,  then  some  more  complex  mov¬ 
ing  wave  is  created  which  does  not  have  a  time-independent  ampli¬ 
tude. 

Waviness  Model 

In  the  previous  report  [8]  it  was  noted  that  the  analytical 
model  for  predicting  performance  of  the  wavy  seal  was  based  on 
assuming  all  negative  fluid  pressures  equal  to  zero.  This  was 
done  in  the  interest  of  simplicity  because  the  problem  solution 
becomes  quite  difficult  for  the  wavy  and  tilted  seal  when  deflec¬ 
tion  must  be  considered.  The  cavitation  solution  model  developed 
previously  [5]  could  not  accommodate  the  tilted  wavy  seal  because 
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the  cavity  shape  may  change  character  from  one  solution  to  the 
next.  Using  the  approximation  It  was  recognized  that  some  error 
would  be  Introduced  Into  the  solution  because  the  correct  cavity 
boundary  conditions  were  not  satisfied. 

Since  the  time  the  previous  report  was  written,  a  paper  has 
been  presented  by  Elrod  [12]  which  provides  a  solution  algorithm 
for  the  cavitation  problem  which  Is  much  more  general  than  those 
developed  previously  [5].  It  provides  relatively  fast,  consis¬ 
tent  solutions  to  the  cavity  problem  using  the  proper  cavity 
boundary  conditions.  This  algorithm  was  incorporated  into  the 
seal  model,  and  previously  given  solutions  wers  recalculated. 

The  algorithm  and  its  Incorporation  Into  the  seal  model  will  now 
be  described. 

Considering  Figure  3-2  which  shows  two  waves  as  an  example, 
the  proper  boundary  conditions  for  the  Reynolds  equations  are  as 


follows.  At  the  Inside  and  outside, 

P (ri)  »  pi  ,  (3-13) 

P(r0)  -  P0  •  (3-14) 

Cavitation  is  most  properly  accounted  for  at  the  upstream  bound¬ 
ary  by 

P  *  Pc  ,  (3-15) 

f  *0  .  (3-16) 

At  the  downstream  boundary,  the  proper  boundary  conditions  are 
considered  to  be  those  given  by  Jakobson  and  Floberg  [13]. 

P  ■  Pc  .  (3-17) 


i£ 

an 


6  cos  a 


& 


{h(2)  ‘  h(D)  • 


(3-18) 
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Rotation  of  Mating  face 


Control  Volume 


Figure  3-2.  Wavy  Seal  Geometry,  Cavitation,  and  Control  Volume 
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where  0  represents  the  downstream  cavity  boundary  and  0 
represents  the  upstream  cavity  boundary.  The  use  of  (3-15), 
(3-16),  (3-17),  and  (3-18)  provides  continuity  across  the 
boundary. 

The  above  equations  were  used  as  stated  in  the  previous  cav¬ 
ity  solutions  [5].  Their  use  requires  that  the  general  shape  of 
the  cavity  be  known  or  the  method  can  be  defeated.  Elrod  [12] 
uses  a  different  approach.  A  dimensionless  variable  a  (e  in  the 
original  paper)  is  introduced  to  provide  continuity.  When  a  is 
less  than  unity,  it  represents  the  fraction  of  flow  in  the  cav¬ 
ity.  When  x  is  greater  than  unity,  it  represents  the  compression 
of  the  fluid.  Thus,  the  same  variable  is  used  throughout  the 
region  and  can  be  solved  for  given  certain  assumptions. 

In  particular,  Elrod  [12]  provides  the  following  relation¬ 
ships: 


(mv) 


x 'convection 


x-lh-l(1 


9-1 )  +  9-^-1 


-1 


(3-19) 


(mx  ^ pressure 
where 


/h3  \ 
\T?n/avg 


Bp, 


9_i(*_i  “  1)  ■  9q  *  (*o  “  1) 


ax 


,  (3-20) 


p  ■  pc  +  s(x  -  1)  and  g  *  1  for  x  _>  1  ,  (3-21) 

p  *  Pc  and  g  *  0  for  x  <  1  .  (3-22) 


S  represents  compressibility,  and  for  this  essentially  incompres¬ 
sible  problem,  a  convenient  arbitrary  value  can  be  used.  The 
subscripts  -1  and  0  represent  node  locations  in  the  x  (e  in  the 
problem  of  interest)  direction  (see  Figure  3-2).  Equation  (3-19) 
applies  only  to  the  tangential  direction.  An  equation  similar  to 
equation  (3-20)  can  be  written  for  the  y  (radial  in  the  problem 
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of  interest)  direction.  Both  equations  (3-19)  and  (3-20)  are 
written  to  describe  flow  across  the  left-hand  boundary  of  a 
control  volume.  Similar  expressions  with  the  subscripts  changed 
apply  to  the  right-hand  boundary. 

Considering  a  complete  control  volume  in  polar  coordinates, 
as  shown  in  Figure  3-2,  the  above  and  similar  expressions  can  be 
used  to  write  a  difference  equation  applicable  to  all  node 
points. 


(Note:  Eauation  (3-23)  is  continued  on  the  next  page.) 
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Equation  (3-23)  represents  the  difference  equation  in  nondimen- 
sional  terms.  The  g-j j  are  as  defined  previously,  "g  is  dimen¬ 
sionless  compressibility.  For  a  grid  of  nodes  ij  (Figure  3-3), 
application  of  equation  (3-23)  results  in  a  system  of  (IMAX-1)* 
(JMAX-2)  linear  equations  in  the  unknown  flow  factors  X-jj. 

Values  of  gij  must  be  assumed  at  the  start.  The  five  diagonal 
coefficient  matrix  of  x-jj  has  a  bandwidth  of  2( JMAX-2)  +  1 
assuming  equations  are  numbered  in  the  r  direction  skipping  from 
end  to  end  to  reduce  bandwidth.  The  matrix  is  nonsymmetrical 
because  of  convective  flow  equation  (3-19).  U  (r^  in  the  problem 
here)  must  be  a  positive  number  or  the  equations  are  invalid. 
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Solution  of  equations  (3-23)  for  the  wavy  seal  was  carried  out 
as  follows. 

1)  All  g-jj  were  set  equal  to  1  at  the  start. 

2)  Values  of  x-jj  on  the  boundaries  were  calculated  using 
equation  (3-21)  and  the  known  pressures. 

3)  The  coefficients  in  equation  (3-23)  were  solved  for  and 
stored  in  a  compacted  form  (Bandwidth  x  number  of  equa¬ 
tions  matrix). 

4)  A  gauss  elimination  and  back  substitution  technique  were 
used  to  solve  for  x-jj. 

5)  Using  (3-21)  and  (3-22),  Pjj  and  g^j  were  solved  for. 

6)  Using  the  p^j,  load  support  and  leakage  were  calculated. 

7)  If  no  convergence,  the  new  g^j  were  used  and  steps  3) 
through  6)  were  repeated. 

Convergence  was  established  by  noting  the  fact  that  the  solution 
eventually  reached  a  stable  oscillation.  That  is,  the  set  of 
values  gij  found  after  a  particular  iteration  would  after  a  few 
iterations  repeat  Itself,  which  means  that  the  solution  cycle 
would  start  over  again.  Once  such  a  repetition  began,  it  was 
concluded  that  the  solution  converged  (to  the  extent  possible). 
Values  of  Pij  used  for  subsequent  calculations  were  based  on 
average  values  within  a  cycle.  It  was  found  that  in  some  cases 
the  g-jj  would  repeat  each  iteration  and  in  other  cases  several 
iterations  occurred  before  the  cycle  repeated.  The  change  in 
load  support  within  a  cycle  was  generally  small.  It  was  con¬ 
cluded  that  the  method  could  not  generally  lead  to  an  absolute 
convergence  because  of  the  finite  grid. 

Using  the  above  technique  a  check  was  made  with  the  solution 
of  the  Findlay  problem  discussed  in  detail  and  solved  in  refer¬ 
ence  [5].  Agreement  with  the  previous  work  was  found  to  be 
excellent,  thus  establishing  complete  confidence  in  the  method. 

The  above  method  does  use  relatively  large  amounts  of  com¬ 
puter  time  compared  to  the  approach  of  using  the  incorrect  cavity 
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boundary  conditions.  However,  it  is  still  fast  enough  that  solu¬ 
tions  to  pressure-caused  deflection  problems  could  be  obtained  at 
a  grid  spacing  of  JMAX=12  and  IMAX=14,  somewhat  coarser  than 
ideal  but  fine  enough  to  obtain  approximate  solutions.  A  typical 
solution  using  the  algorithm  was  obtained  for  the  base  case  prob¬ 
lem  in  this  work  with  an  ideal  combination  of  wav i ness  and  tilt 
at  n  *  9.  The  case  presented  is: 

Pi  3  0 

p0  =  3.45  MPa 
Pm  3  262  MPa 
Ps  3  0*1  tyn 
psp  *  0.207  MPa 

C  3  0.51  urn  (3-24) 

B  =  1.0 

n  «  6.83  •  10-4  Pa  •  s 
rpm  *  1800 
r -j  =  48.26  rrm 
r0  =  54.04  rim 
rc  3  49.05  rim 
n  3  9 

Nondimensional  constants  are: 

Pq  3  2.46  •  10"3 
Pm  3  1-87  •  10-1 
ps  3  1.87  •  10-2 

(3-25) 

psp  3  1.48  •  10-4 
rQ/c  3  1.044  .  105 
3  0.910 
7C  =  0.925 
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Wav i ness  selected  is: 


v0  a  0.36 
To  *  24 


(3-26) 


Using  the  new  algorithm,  the  computed  values  are: 


1800  rpm 
u  -  0.0019 
Q  »  2.3  ml/min 
T  »  0.54  N  •  m 
X  *  100 
hQ  *  1.23 


1  rpm 
u  -  0.0277 
Q  »  0.6  ml/min 
T  -  7.8  N  •  m 
%  »  72.0 
hQ  =  0.66 


(3-27) 


Figure  3-4  shows  the  film  thickness  distribution  and  Figure  3-5 
shows  the  pressure  distribution  for  1800  rpm.  The  cavity  is 
identified  as  the  region  of  zero  pressure. 

The  new  algorithm  has  been  used  for  calculating  all  of  the 
theoretical  results  shown  in  Chapter  4. 


Optimum  Seal  Design  and  Number  of  Waves 

Since  the  previous  report  was  written  [8],  further  work  on 
optimum  seal  design  and  number  of  waves  has  been  done  and  pre¬ 
sented  in  reference  [14].  This  work  was  carried  out  using  the 
new  algorithm  above.  The  essential  features  of  this  work  will  be 
given  here. 

Considering  the  base  case  values  of  (3-24)  and  optimum  con¬ 
ditions,  reference  is  made  to  Figure  3-7  where  dynamic  perform¬ 
ance  is  shown  as  a  function  of  To  f°r  9i'ven  value  of  70.  Con¬ 
sider  first  the  case  for  70  =  5.0  which  is  also  shown  by  Figure 
3-6.  As  To  increases  leakage  begins  to  decrease  rapidly  and 
torque  slowly  begins  to  increase.  In  the  static  case  shown  in 
Figure  3-8,  there  is  a  distinct  minimum  in  the  torque  curve  at  a 
particular  value  of  <j>.  This  value  corresponds  to  the  middle 
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Figure  3.5.  Pressure  Distribution  -  p  x  10" 


Figure  3-6.  Worn  In  Film  Shape  at  Various  t  With  fi 


film  thickness  configuration  shown  by  Figure  3-6,  which  is  the 
condition  where  touching  just  begins  to  occur  completely  around 
the  seal  at  rj.  As  discussed  previously  [8],  this  condition 
can  be  shown  to  be  where 

1  •  -  V°  —  (3-28) 

°opt1mum  rc  *  r1 

For  7C  *  5  and  70  »  0.99,  7o0ptimum  a  62.5.  Understandably  as  7o 
is  increased  beyond  this  point,  leakage  becomes  quite  low  both 
statically  and  dynamically  because  a  continuous  sealing  band  is 
formed.  It  is  also  clear  that  as  "$o  increases  the  friction  torque 
begins  to  increase  for  the  static  case  because,  as  shown  in  Fig¬ 
ure  3-6,  a  portion  of  the  converging  film  which  reduces  hydro¬ 
static  load  support  Is  lost.  For  "$o  values  lower  than  the 
point  above,  leakage  is  high  because  of  large  gaps,  and  static 
torque  is  high  because  the  effect  of  taper  is  lost  at  large  gaps. 

This  illustration  shows  how  a  combination  of  waviness  and 
tilt  leads  to  better  performance  than  waviness  alone.  With  the 
proper  combination,  leakage  can  be  much  smaller  because  there  are 
no  large  gaps  as  caused  by  parallel  face  waviness.  Friction, 
particularly  static  friction,  can  be  made  much  lower  than  for 
radially  parallel  waviness  because  the  tilt  introduces  a  converg¬ 
ing  taper  over  much  of  the  face  of  the  seal. 

For  n  *  9,  70  *  5  is  not  practical  as  will  be  shown  later 
and  has  been  used  to  demonstrate  the  conditions  where  performance 
approaches  an  optimum.  The  case  where  70  *  0.36  illustrated  in 
Figures  3-7  and  3-8  is  physically  possible  but  represents  a  de¬ 
parture  from  the  behavior  of  the  previous  case.  The  optimum 
point  is  at  To  =  24.  In  Figure  3-7  it  is  clear  that  friction 
is  high  for  70  less  than  24.  However,  friction  and  leakage 
both  continue  to  decrease  with  increasing  To  far  beyond  the  opti¬ 
mum  point.  The  reason  for  this  behavior  is  that  the  seal  radius 
at  which  continuous  contact  occurs  (see  Figure  3.6c)  moves  in 
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very  slowly  with  respect  to  <f>  for  this  case  as  compared  to  the 

previous  case.  Thus,  some  slight  improvement  is  found  beyond  the 

optimum  point  whereas  in  the  previous  case  the  change  to  less 

advantageous  operating  conditions  occurs  so  rapidly  the  potential 

slight  improvement  with  $  beyond  the  optimum  point  is  hardly 

noticed.  In  any  case,  the  designer  can  choose  any  physically 

achievable  operating  point  beyond  the  optimum  point  where 

** 

performance  is  best  or  where  performance  is  most  satisfactory. 

The  selection  of  70  depends  to  a  large  extent  on  what  is  physi¬ 
cally  achievable  as  well  as  acceptable  maximum  leakage.  Accord¬ 
ing  to  reference  [8]  for  a  three-wave  study,  for  the  optimum 
point  there  is  a  tradeoff  between  leakage  and  wear  rate  which 
depends  on  70.  To  achieve  maximum  benefit  one  should  choose  an 
operating  point  with  less  than  desired  leakage  and  with  %  fluid 
film  load  support  such  that  the  ratio  0.5/(l  -  F)*  gives  a  suita¬ 
ble  wear  reduction  factor  compared  to  parallel  face  performance. 
It  is  also  desirable  that  static  torque  be  minimized  so  that  low 
speed  operation  or  startup  will  not  damage  the  seal. 

The  means  selected  in  this  work  to  produce  alternating  tilt 
and  waviness  is  to  apply  amoving  cosinusoidal ly-distributed 
moment  to  the  ring,  i.e., 

m9  =  m3o  cos(n9  +  flt)  (N  •  m/m)  ,  (3-29) 


where  the  moment  must  be  expressed  as  moment/unit  length  of  cir¬ 
cumference.  From  reference  [8],  it  can  be  shown  that  the  above 
moment  produces  a  cosinusoidal  v  and  $  with  amplitudes  given  by 


v 


o 


(1  +  A) 


(3-30) 


*The  0.5  refers  to  a  50  percent  fluid  pressure  load  support  for  a 
parallel  face  balance  ratio  one  seal. 
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(3-31) 


1 


%rc  1  +  An2 


0  “x  (n  -  l)' 

The  physical  forces  producing  the  moment  need  not  be  liter¬ 
ally  cosinusoidal ly  distributed.  Any  means  of  applying  concen¬ 
trated  loads  or  discrete  length  distributed  loads  whose  nth  term 
Fourier  series  representation  corresponds  to  equation  (3-29)  is 
sufficient.  For  example,  in  the  device  shown  in  Figure  3-1,  the 
loads  are  derived  from  gas  pressure  acting  over  a  circular  area. 
The  Fourier  analysis  of  the  staggered  pattern  of  the  holes  gives 


sin7FT 


(3-32) 


Higher  harmonic  terms  are  also  present,  but  these  generally  give 
negligibly  small  deflections.  On  the  other  hand,  one  must  be 
certain  that  unwanted  lower  harmonics  are  not  introduced  by  what¬ 
ever  method  is  used  to  apply  force. 

Now  if  the  method  selected  for  applying  the  moment  is  that 
above,  to  obtain  the  optimum  ratio  of  <j>0  to  v0,  substitution 
of  equations  (3-30)  and  (3-31)  into  (3-28)  gives 


1  +  An* 


(3-33) 


For  the  data  given  by  (3-25),  this  would  require  that  A  =  3.14, 
a  realistic  value.  Clearly,  if  A  does  not  satisfy  this  require¬ 
ment,  or  any  other  non-optimum  case  requirement  involving  the 
ratio  of  equations  (3-30)  and  (3-31),  then  the  geometry  must  be 
adjusted,  changing  rc  and  A,  until  equation  (3-28)  is  satis¬ 
fied.  In  general,  the  ratio  $0/v0  has  definite  limits  imposed 
by  the  method  of  application.  Some  design  experimentation  must  be 
performed  both  in  choosing  a  ratio  4>0/v0  which  gives  suitable 
performance  and  which  is  physically  achievable. 
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In  practical  applications,  there  are  errors  of  approximation 
associated  with  the  calculations  of  both  rc  and  A  such  that  a 
less  than  optimum  condition  might  in  fact  result.  To  avoid  this 
possibility,  it  has  been  determined  that  the  contact  point  as 
shown  in  Figure  3-6  should  be  designed  to  be  at  some  radius 
greater  than  the  inside  radius  as  shown  by  3-6c.  The  condition 


-  —  +  r 


(3-34) 


defines  the  contact  radius.  If  the  design  contact  radius  rd  is 
specified  to  be  some  fraction  f  of  the  face  width  outside  of  the 
inside  radius,  then 


=  r  •  +  f  ( r  -  r  • ) 

i  '  o  i ; 


(3-35) 


Then  the  required  ratio  for  v0/$Q  is 


r  =  rc  ■  ri  -  ,(ro  -  ri> 


(3-36) 


and  equations  (3-30)  and  (3-31)  may  be  applied  to  determine  the 
necessary  geometry  to  provide  for  this  more  realistic  design 
condit  ion . 

The  number  of  waves  to  be  used  is  not  arbitrary,  particu¬ 
larly  if  the  seal  is  to  operate  under  various  operating  condi¬ 
tions.  It  is  shown  in  reference  [8]  that  the  deflection  due  to 
face  pressure,  referred  to  above,  changes  with  operating  pressure 
and  speed.  This  means  that  the  faces  must  wear  into  a  new  equi¬ 
librium  shape  for  each  change  of  conditions.  This  will  acceler¬ 
ate  wear  and  is  undesirable. 

To  counter  this  problem,  it  is  possible  to  design  the  seal 
ring  so  that  the  desired  deflection  is  large  relative  to  the 
changing  face  pressure  deflection.  Thus,  when  operating  condi¬ 
tions  change,  there  will  be  no  significant  change  in  the 
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equilibrium  face  profile.  This  objective  can  be  achieved  by  not¬ 
ing  the  results  of  the  ratios  of  the  face  pressure-caused  deflec¬ 
tion  to  the  desired  deflection  as  given  by  equations  (4-47)  and 
(4-68)  of  reference  [8]. 
v. 

P.~.  .  .  e. 

A  +  n  (3-37) 


o  r 
^o__c 


om 


m 


9 


90 


n‘(A  +  1) 


?opt 


rom. 


o  r 
Hyo  c 


m 


90 


A  +  1 
An2  +  1 


(3-38) 


As  n  becomes  large,  the  second  of  these  two  ratios  becomes  small. 
This  means  that  the  tilt  angle  40  is  not  significantly  Influ¬ 
enced  by  face  pressure  py  and  therefore  not  Influenced  signifi¬ 
cantly  by  changing  operating  conditions.  v0  is  not  greatly 
affected  by  n.  However,  it  is  the  change  In  with  operating 
conditions  that  leads  to  accelerated  wear.  Thus,  by  using  large 
n,  the  detrimental  effect  of  unwanted  deflection  can  be  minimized. 
For  the  particular  example  under  study  with  n  »  9,  changes  in  <f> 
due  to  changes  in  operating  conditions  were  less  than  one  percent. 

The  above  considerations  have  been  made  for  the  design  of  a 
new  experimental  seal  using  n  a  9  and  the  off-optimum  condition 
mentioned.  It  is  expected  that  the  new  design  will  be  insensitive 
from  a  wear  standpoint  to  the  operating  pressure  and  speed. 


Offset -Land  Theory 

The  offset-land  seal  is  shown  in  Figure  3-9.  The  offset 
portion  of  the  seal  Is  continuously  machined  away  by  protrusions 
attached  to  the  hard  face  as  described  in  the  previous  chapter. 

To  calculate  performance  of  the  offset- land  seal,  a  computer  pro¬ 
gram  developed  previously  [9]  was  modified  to  accommodate  the  off¬ 
set  geometry.  Leakage,  load  support,  and  friction  torque  are  cal¬ 
culated  as  usual.  To  compare  the  anticipated  wear  rate  of  the 
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offset  design  to  a  conventional  parallel  face  design,  the  mechan¬ 
ical  contact  pressure,  to  which  wear  rate  is  proportional,  was 
calculated.  Thus,  the  relative  wear  rate  is  given  by 


RWR  * 


iWoffse‘> 


(3-39) 


When  this  ratio  is  1.0,  the  wear  ratio  of  the  land  of  the  offset 
land  is  the  same  as  that  of  a  parallel  face  seal.  As  the  ratio 
becomes  smaller,  the  wear  rate  of  the  offset -land  seal  would  be 
expected  to  be  the  RWR  fraction  of  the  parallel  face  wear  rate. 

Performance  was  evaluated  over  a  wide  range  of  variables. 

It  was  determined  that  the  concept  is  useful  only  for  seals  hav¬ 
ing  a  balance  ratio  less  than  1.0.  The  following  studies  for  B  * 
0.75  were  based  on  the  following  parameter  values. 

Pi  *  0 

p0  *  3.45 .MPa 
Pn,  *  262  MPa 
Ps  3  0-1  Pm 
psp  =  0.207  MPa 

a  =  0.17  um  (3-40) 

n  =  6.83  •  lO’4  Pa  •  s 
rpm  =  1800 
r^  *  48.26  mm 
r0  =  54.04  mm 
^offset  -  variable 
A/B  -  variable 


Also, 

^offset  =  hoffset^a  '  (3-41) 


Figure  3-10  shows  the  effect  of  land  width  on  performance  holding 
^offset  3  5  (h0ffSet  3  0*85  w).  From  the  standpoint  of  minimiz¬ 
ing  both  torque  and  relative  wear  rate,  there  appears  to  be  an 
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Figure  3-10.  Effect  of  Land  Width  on  Performance. 


optimum  land  ratio  around  0.2.  Wear  rate  becomes  less  than  0.1 
that  of  a  comparable  parallel  face  seal.  Leakage  values  at 
^offset  “  5  are  nodest. 

Figure  3-11  shows  the  effect  of  h0ffset  on  performance 
with  land  ratio  held  constant  at  0.2.  The  relative  wear  rate 
quickly  drops  to  insignificant  values  as  ^0ffset  increases  beyond 
8.  Further  increases  in  h0ffset  do  little  to  decrease  torque 
and  cause  leakage  to  increase  to  large  values.  Based  on  these 
studies,  a  good  design  point  would  be  for  an  ^Qffset  between  6 
and  10  (1.0  -  1.7  um).  This  would  result  in  wear  rates  at  least 
a  factor  of  10  lower  than  parallel  face  seals  and  low  torque 
while  leakage  would  be  less  than  10  ml/min. 

The  results  of  the  theory  described  here  are  compared  to 
experimental  results  in  the  next  chapter. 
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Chapter  4 

Comparison  of  Theoretical  and  Experimental  Results 
Three-Wave  Seal  Tests 

Using  the  new  algorithm  described  in  the  previous  chapter, 
the  theoretical  results  for  the  base  case,  three-wave,  moving 
waviness  test  as  a  function  of  waviness  pressure  were  recomputed. 
These  results  are  shown  in  Table  4-1.  The  results  are  similar  to 
those  in  Table  4-4  in  the  previous  report  [8]  which  were  obtained 
using  an  approximate  method.  The  two  columns  for  hpj,  and  hpbcaic 
have  signs  opposite  to  those  in  the  previous  results.  This 
occurs  because  the  direction  of  rotation  was  reversed  for  the 
present  calculations  as  is  required  by  the  new  algorithm.  The 
other  significant  difference  is  that  the  predicted  leakage  values 
are  lower  than  the  previous  predictions.  The  leakage  is  compared 
to  experimental  results  in  Figure  4-1.  Agreement  is  better  than 
that  obtained  using  the  previously  calculated  results.  Torque 
values  were  adjusted  by  reducing  the  friction  coefficient  from 
the  value  used  for  computation.  These  values  are  compared  to 
experimental  results  in  Figure  4-2.  Agreement  is  reasonable  and 
similar  to  that  obtained  previously. 

In  summary,  the  use  of  the  more  correct  algorithm  makes  some 
improvement  in  predicting  leakage  for  the  three-wave  seal.  Al¬ 
though  the  difference  here  is  not  large,  it  was  found  for  the 
nine-wave  case  that  the  approximate  model  was  totally  in  error 
because  of  the  extensive  cavitation  in  that  case.  Without  the 
use  of  the  new  algorithm,  total  disagreement  for  the  nine-wave 
case  was  found. 

Radial  Taper 

Using  the  theoretical  model  presented  in  detail  in  reference 
[9],  calculations  were  made  to  predict  performance  for  the  radial 
taper  tests.  The  model  predicts  both  short-term  (initial)  and 
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long-term  performance.  Table  4-2  presents  typical  computed  re¬ 
sults  for  B  =  0.75  and  a  roughness  of  a  =  0.54  un  which  was  the 
average  final  roughness  for  these  tests.  The  computed  values 
show  a  reduction  in  leakage  and  an  increase  in  torque  as  the 
radial  taper  goes  from  convergent  to  divergent  as  previously 
discussed. 

Considering  first  the  initial  performance.  Figures  4-3  and 
4-4  show  predicted  torque  and  leakage  as  functions  of  radial 
taper  for  the  B  =  0.75  cases  compared  to  experimental  results. 

Two  theoretical  curves  are  given,  one  for  a  =  0.54  un  as  dis¬ 
cussed  above  and  one  for  a  =  0.15  urn.  The  reason  for  the  second 
curve  will  become  apparent  later.  The  data  on  both  figures 
agrees  reasonably  well  with  theory  at  a  -  0.54  thus  substantiat¬ 
ing  the  model. 

Figures  4-5  and  4-6  show  a  similar  comparison  made  for  the 
B  =  1.0  tests  whose  results  were  reported  in  the  previous  report 
[8].  Figure  4-5  shows  good  agreement  between  theory  and  experi¬ 
ment  at  a  =  0.54.  Figure  4-6  shows  that  a  =  0.54  greatly  over¬ 
estimates  leakage,  a  =  0.15  provides  a  better  appreciation  for 
leakage  (which  was  unmeasurable  on  three  of  the  four  tests). 

This  discrepancy  points  out  the  reason  for  advancing  the 
socket  theory  in  Chapter  2.  It  is  clear  from  Figure  4-6  (and 
other  data  to  be  discussed)  that  using  the  actual  roughness  can 
overestimate  the  leakage  in  many  cases.  The  effective  value  of 
the  roughness  may  be  much  smaller  than  the  actual  roughness  as 
explained  previously.  What  is  not  clear  is  why  the  higher  value 
of  roughness  more  accurately  predicted  the  results  in  the  B  = 

0.75  cases.  It  is  clear  that  the  theory  is  inadequate.  Neither 
the  original  theory  nor  the  pocket  theory  are  consistent.  One 
thing  that  is  known  is  that  the  surface  roughness  changes  during 
the  course  of  a  test.  This  may  in  part  explain  these  differences 
in  that  the  surfaces  evolve  differently  in  the  two  types  of  tests 
because  the  mechanical  contact  pressures  are  different.  While 
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Table  4-2.  Theoretical  Results— Radial  Taper 
a  =  0.54  un  (21.25  yin.),  B  »  0.75 


<t>i 

( um/m) 

Initial  Torque 
N*m  (in. -lb) 

Final 

N»m 

Torque* 

( in. -lb) 

Leakage 

9  1  hr 
(cm3/min) 

Final 

Leakage* 

(cm3/min) 

1013 

0.2 

(1.5) 

0.19 

(1.7) 

48.5 

41.4 

613 

0.8 

(6.6) 

8.12 

(71.9) 

18.0 

4.18 

40 

3.5 

(31.3) 

8.18 

(71.4) 

6.43 

4.17 

-320 

6.9 

(60.6) 

8.17 

(72.3) 

4.25 

4.17 

-750 

11.5 

(101.7) 

8.18 

(72.4) 

3.75 

4.17 

-1110 

15.0 

(132.4) 

8.18 

(72.4) 

3.79 

4.17 

♦Based  on  length  of  corresponding  test. 
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the  theory  does  predict  leakage  trends,  much  more  work  will  be 
needed  to  understand  the  nature  of  surface  roughness  as  it 
affects  leakage  in  order  to  accurately  predict  leakage  quantita¬ 
tively. 

Figures  4-7  and  4-8  show  torque  and  leakage  versus  time  for 
a  highly  convergent  B  =  0.75  case.  The  experimental  data  are 
compared  to  theoretical  data  computed  based  on  the  time-dependent 
mode  of  the  radial  taper  model  [9].  The  experimental  results 
are  scattered  and  there  is  no  discernable  trend.  The  theoretical 
values  change  very  little  over  this  time  interval,  and  one  can 
say  that  theory  and  experiment  agree  to  some  extent. 

In  Figures  4-9  and  4-10,  the  inadequacy  of  the  time- 
dependent  model  becomes  apparent  for  this  slightly  less  converg¬ 
ing  case.  In  Figure  4-10  for  o  =  0.54  urn,  it  is  predicted  that 
some  of  the  converging  taper  will  wear  away  causing  the  torque  to 
increase  and  leakage  to  decrease  with  time.  This  is  clearly  not 
the  case.  Using  o  =  0.15  urn  helps  very  little.  Torque  agrees  in 
that  there  is  virtually  no  change  with  time,  but  predicted  leak¬ 
age  is  now  too  low.  It  is  likely  that  not  only  does  the  model 
suffer  from  an  inadequate  characterization  of  the  surface  rough¬ 
ness  as  discussed  previously,  it  also  suffers  because  the  wear 
model  used  is  not  sufficiently  accurate,  even  though  the  wear 
constant  used  is  based  on  experimental  data.  In  general,  it  has 
been  found  for  the  converging  cases  that  very  little  wear  occurs 
during  a  test,  much  less  than  predicted  based  on  the  model  and 
parallel  face  test  results. 

Figures  4-11  and  4-12  are  for  a  highly  divergent  case.  Here 
in  Figure  4-11  some  agreement  is  found  in  that  experimental  tor¬ 
que  does  decrease  with  time,  although  not  as  much  as  predicted. 
Measured  leakage  again  agrees  much  better  with  the  a  =  0.15  urn 
analysis.  Figure  4-12  shows  that  in  the  long  run  torque  varies 
widely  with  time  which  is  characteristic  of  a  parallel  face  test. 
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©  Experimental  Results 
TIME  (h)  —  Theory,  a  ■  .15  ym 

Theory,  a  =  .54  ym 

Figure  4-7.  Torque  and  Leakage  Versus  Time,  Test  64  -  B  =  0.75. 
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Figure  4-9.  Torque  and  Leakage  Versus  Time,  Test  65  -  B  =  0.75. 
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Figures  4-13  and  4-14  for  two  tests  at  B  =  1.00  illustrate 
further  the  problems  with  the  model.  Wear  is  predicted  to  occur 
much  faster  than  it  actually  does.  The  wear  constant  used  is 
clearly  not  a  constant  but  apparently  is  a  strong  function  of  the 
mechanical  contact  pressure  and  other  factors. 

When  operation  is  such  that  parallel  face  operation  results 
over  a  sufficient  portion  of  the  face  that  a  worn  in  slope  can  be 
determined,  then  one  may  determine  the  radial  taper  or  rotation 
caused  by  thermal  effects.  Figures  4-15  and  4-16  show  experimen¬ 
tally  determined  radial  taper  compared  to  predicted  thermal  taper 
as  a  function  of  final  measured  torque.  These  results  were  not 
based  strictly  on  the  radial  taper  tests  but  included  a  number  of 
originally  parallel  face  cases  particularly  for  the  B  *  1.00 
cases.  The  lower  experimental  values  for  B  =  1.00  correspond 
well  to  the  theory.  The  B  a  0.75  values  indicate  that  the  model 
underestimates  the  thermal  taper.  This  is  thought  to  be  a  result 
of  the  fact  that  the  model  uses  film  coefficients  based  on 
experiments  which  are  for  a  geometry  somewhat  different  than  that 
of  a  seal.  By  plotting  taper  as  a  function  of  torque,  errors 
caused  by  imprecise  friction  coefficient  estimates  are  elimina¬ 
ted.  In  spite  of  the  difference,  the  model  does  provide  a  good 
first  approximation  to  predicting  thermal  rotation. 

High  Temperature  Tests 

As  discussed  previously,  vrfien  a  seal  is  operated  at  a  tem¬ 
perature  such  that  the  seal  face  temperature  is  above  the  satura¬ 
tion  temperature  of  the  fluid  at  the  discharge  pressure,  perform¬ 
ance  may  change  dramatically.  In  previous  work  [11],  a  model  was 
developed  to  predict  this  behavior.  The  model  has  been  applied 
to  the  test  seal  geometry. 

The  details  of  the  geometry  are  shown  in  Figure  4-17.  The 
parameter  values  for  the  computer  solutions  are  as  follows. 
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Figure  4-13.  iorque  and  Leakage  Versus  Time,  Test  24 


Seal  Geometry  for  Two-Phase  Model . 


Seal  Parameters: 


ri 

rl 

r2 

r3 

n 

s2 

s3 

«4 

**4 

r5 

s5 

Psp 

Pm 

Ps 


53.04 

nm 

(2.088 

in.) 

48.26 

nm 

(1.900 

in. ) 

54.28 

mm 

(2.137 

in.) 

59.69 

ran 

(2.350 

in.) 

48.26 

mm 

(1.900 

in. ) 

for  B  *  1.0; 

49.49 

nm 

(1.949 

in.) 

for  B  »  0.75 

3.18  nm  i 

(0.125  ' 

in.) 

13.34 

mm 

(0.525 

in.) 

15.24 

mm 

(0.600 

in.) 

0.0  mm  for  B  *  1.0; 

6.45  rrm  (0.254  in.)  for  B  *  0.75 
47.19  mu  (1.858  in.) 

57.14  mn  (2.250  in.) 

19.96  mm  (0.786  in.) 

2.07  x  10s  N/m2  (30.0  lb/in. 2) 

2.62  x  10®  N/m2  (38000  1b/in.2) 
2.62  x  107  N/m2  (3800  lb/in. 2) 
0.5  tfti  (19.69  yin.) 


Heat  Transfer  Parameters: 

H0  *  (given  later) 

H-,-  *  (given  later) 

kj  =  11.25  W/m*K  (6.5  Btu/hr-ft-’F) 

k2  *  70.00  W/m*K  (40.0  Btu/hr-ft-'F) 

Operating  Conditions: 

p0  *  1.72  x  106  N/m2  (250  lb/in.2  absolute) 

Pi  *  8.43  x  10^  N/m2  (local  atmospheric  pressure 

(12.2  Ib/in.2) 
to  *  188  1/s  (1800  rev. /min) 


Fluid  Properties: 

RH2o  *  460  N«m/kg*K  (0.110  Btu/lb^.’F) 

Vv  Mg,  Pv  Psat,  Hf ,  Hfg  —  from  steam  table  data. 
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The  value  selected  for  ps  gives  the  mechanical  friction 
coefficient  F  =  0.1  when  all  of  the  load  is  supported  by  mechani¬ 
cal  contact.  Seal  ring  1  properties  of  k^  and  pm  are  typical  of 
a  carbon  material.  The  thermal  conductivity  of  seal  ring  2,  k2» 
is  typical  of  tungsten-carbide.  The  RMS  combined  roughness  a  is 
the  approximate  roughness  value  that  evolves  during  test  opera¬ 
tion. 

Two  different  convection  coefficients  HQ  and  have  been 
included  in  the  heat  transfer  model.  H0  is  approximated  by 
using  the  convection  coefficient  of  a  horizontal  cylinder  rotat¬ 
ing  in  a  tank  of  water  [15].  By  using  the  formula 

NN(j  =  0.133  x  Re2/3  x  Pr1/3  (4-1) 

the  convection  coefficient,  H0,  for  different  operating  condi¬ 
tions  can  be  determined.  The  H„  values  used  for  prt  *  1.72  x 

6  2  0  0 
10  N/m  and  different  operating  temperatures  are  shown  in  Table 

4-3.  Oue  to  the  complexity  of  the  seal  assembly,  it  is  not  easy 
to  determine  the  H-j  value  accurately.  An  approximate  value  of 
Hi  =  7000  w/m  K  is  used  in  computation.  As  discussed  later, 
this  value  is  justified  by  investigating  the  equivalent  convec¬ 
tion  coefficients  on  different  seal  parts  and  then  making  an 
overall  approximation  for  Hi. 

Using  these  data  and  the  model,  theoretical  performance  for 
test  seal  torque  has  been  calculated  and  is  compared  to  experi¬ 
mental  results  in  Figures  4-18  and  4-19  for  B  *  1.00  and  B  « 

0.75.  Two  values  of  mechanical  friction  coefficient  were  used 
for  the  calculation  as  shown.  At  F  =  0.15,  agreement  with  theory 
is  reasonable  up  to  the  transition  point.  The  approximate  tem¬ 
perature  of  the  transition  point  is  predicted.  Beyond  the  tran¬ 
sition  point,  experimental  friction  torque  is  much  higher  than 
predicted.  It  is  thought  that  this  is  probably  caused  by  ex¬ 
tremely  high  coefficients  of  mechanical  friction  associated  with 
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Table  4-3.  Outside  Convection  Coefficients,  Hg 
p0  =  1.72  MPa 


T, 

C  c) 

CF) 

H0(Btu/hr.ft2°F) 

H0( W/m2*K) 

37.8 

100 

3120 

17700 

65.6 

150 

3980 

22600 

93.3 

200 

4400 

25000 

121.1 

250 

4670 

26500 

148.9 

300 

5070 

28800 

176.7 

350 

5350 

30400 

182.2 

360 

5370 

30500 

187.8 

370 

5400 

30700 

193.3 

380 

5420 

30800 

198.9 

390 

5440 

30900 
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Figure  4-13.  Torque  as  a  Function  of  Temperature  -  B  *  1 .0 
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Figure  4-19.  Torque  as  a  Function  of  Temperature, 
3  =  0.75. 
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operation  of  the  faces  in  essentially  superheated  vapor.  Exper¬ 
iments  are  needed  to  verify  this  hypothesis;  such  data  is  not 
available  in  the  literature. 

It  can  be  shown  theoretically  that  the  temperature  differ¬ 
ence,  *  aT,  is  proportional  to  the  friction  torque  for 

different  operating  temperatures.  As  shown  in  Figure  4-20, 
experimental  results  also  show  this  trend,  although  theory  over¬ 
estimates  the  aT.  Thus,  the  AT  data  can  be  used  as  an  indication 
of  friction.  In  Figures  4-21  and  4-22,  such  experimental  results 
are  compared  to  theory.  The  experimental  values  of  aT  are  aver¬ 
aged  over  the  last  few  hours  of  the  test.  The  theoretical  tem¬ 
perature  estimate  is  based  on  the  calculated  temperatures  at  the 
point  in  the  seal  body  where  temperature  is  measured. 

Considering  Figure  4-18  and  Figure  4-21  (or  Figure  4-19  and 
Figure  4-22),  the  prediction  of  the  effect  of  operating  tempera¬ 
ture  on  friction  is  consistent  between  that  based  on  aT  and  that 
based  on  torque.  Given  that  aT  is  a  reliable  measurement,  this 
also  suggests  that  the  measurement  of  torque  is  satisfactory. 

Experimental  results  in  Chapter  2  show  that  after  a  running 
in  period  no  leakage  has  ever  been  measured  for  the  high  tempera¬ 
ture  tests.  Based  upon  the  theoretical  analysis,  the  leakage 
rate  is  determined  primarily  by  the  surface  roughness  (a  cubic 
relationship).  The  standard  deviation  of  surface  roughness  of 
the  carbon  after  test  was  found  to  vary  from  0.4  to  1.0  um. 
Leakage  results  shown  in  Figure  4-23  are  based  on  a  -  0.5  um, 
which  compares  favorably  to  measured  values.  Since  no  leakage 
has  been  measured,  it  is  thought  that  the  effective  roughness  is 
lower  than  the  measured  value  as  discussed  with  regard  to  the 
radial  taper  model.  The  pocket  theory  previously  discussed  sug¬ 
gests  that  flow  is  controlled  primarily  by  the  roughness  magni¬ 
tude  associated  with  the  flat  peaks.  This  value  amounts  to 
approximately  0.1  um  of  standard  deviation,  much  less  than  the 
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Figure  4-20.  aT  Versus  Torque,  B  *  1.0,  F  Varies 
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aT  as  a  Function  of  Temperature  -  S  = 
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Figure  4-22.  AT  as  a  Function  of  Temperature,  B  =  0.75. 
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standard  deviation  measured  by  taking  the  entire  surface  into 
account.  By  using  the  value  a  =  0.1  um,  the  computed  results  of 
leakage  are  less  than  0.05  cm^/min.  At  this  low  rate  it  is 
expected  that  water  will  evaporate  from  the  leakage  counter 
before  it  can  be  registered.  Thus,  the  experimental  results  are 
consistent  with  the  theory  as  modified. 

Considering  Figures  4-18  and  4-19,  it  has  been  shown  that 
the  higher  mechanical  friction  coefficient  (F  =  0.15)  may  be  used 
such  that  the  theoretical  curve  of  torque  intersects  the  experi¬ 
mental  data.  In  Figure  4-21,  it  can  be  seen  that  the  theoretical 
aT  value  based  on  this  mechanical  friction  coefficient  is  some¬ 
what  higher  than  the  experimental  results.  This  indicates  that 
the  validity  of  the  heat  transfer  model  must  be  examined  more 
closely. 

Figure  4-24  shows  the  geometry  of  the  actual  seal  assembly 
and  that  of  the  heat  transfer  model .  It  can  be  seen  that  the 
geometry  of  the  actual  seal  is  somewhat  different  from  that  of 
the  heat  transfer  model  particularly  in  region  D.  However,  modi¬ 
fication  of  the  geometry  of  the  heat  transfer  model  to  fit  the 
actual  seal  geometry  has  been  shown  to  be  unnecessary  since  the 
theoretical  temperature  distribution  in  the  portion  0  is  very 
close  to  T*. 

In  the  following  paragraphs,  the  selection  of  heat  transfer 
parameter  values  for  the  B  =  1.0  seal  and  the  sensitivity  of  Tf 
to  the  change  of  each  parameter  value  are  examined. 

1.  During  operation,  the  W-C  ring  is  rotating  while  the 
carbon  ring  remains  stationary.  The  calculated  value  of 
H0  using  equation  (4-1)  is  a  convection  coefficient  of 
a  rotating  cylinder  in  water.  Thus,  there  might  be  an 
error  in  the  convection  coefficient  for  the  stationary 
carbon  ring.  However,  considering  the  flow  around  the 
seal  rings  and  the  complexity  of  the  seal  assembly,  the 
assumption  used  might  be  the  best  approximation  which 
can  be  made. 


137 


\ 

\ 

l 


Heat  Transfer  Model 


Figure  4-24.  Geometry  of  Actual  Seal  Assembly  and  Heat 
Transfer  Model . 


2.  The  thermal  conductivities,  k^  and  kg,  are  typical  val¬ 
ues  for  carbon  and  W-C.  However,  they  are  not  actually 
measured  values.  It  is  known  that  the  value  of  kg  is  in 
the  range  of  43  to  86  W/m*K.  Based  upon  theoretical 
analysis  using  kg  *  86  W/m*K,  the  aT  curve  is  generated 
as  shown  in  Figure  4-25  as  curve  c.  Curve  b  is  based  on 
the  value  k  =  70  W/m*K  assumed  for  the  previously  given 
results.  The  result  shows  that  the  aT  prediction  is 
quite  sensitive  to  the  value  of  k£. 

3.  Due  to  the  complexity  of  seal  assembly,  it  is  not  easy 
to  accurately  determine  the  inside  convection  .  Con¬ 
sidering  Figure  4-24,  it  is  believed  that  positions  1, 

3,  and  4  are  filled  with  water  with  temperature  Tm  dur¬ 
ing  operation  when  operation  is  below  the  saturation 
temperature.  By  investigating  the  thermal  conductivity 
of  water  and  the  gap  thickness  between  the  seal  and  the 
surrounding  components  and  assuming  these  components  are 
at  temperature  Tw,  equivalent  (thermal  conductivity 
of  water  divided  by  the  gap  thickness)  can  be  deter¬ 
mined.  The  Hi  value  at  position  2,  assuming  air  is 
around  this  area,  is  determined  by  using  equation  (4-1). 
After  H.j  values  on  different  seal  parts  have  been 
obtained,  an  overall  approximated  value  of  H_-  .=  7000 
W/m K  was  selected  for  computation.  In  order  to 

observe  the  effect  of  the  value  on  the  aT  predic- 

1  2 

tion,  a  aT  curve  based  on  H.  s  9000  W/m K  is  shown 
in  Figure  4-25  as  curve  d.  Results  show  that  theore¬ 
tical  prediction  of  aT  is  quite  sensitive  to  the  Hi 
value. 

Thus,  for  more  accurate  temperature -prediction,  further  in¬ 
vestigation  of  these  heat  transfer  parameters  is  required.  How¬ 
ever,  considering  the  numerous  approximations  and  uncertainties  in 
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Figure  4-25.  Heat  Transfer  Parameters  -  B  =  1.0. 
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both  theory  and  experiment,  the  heat  transfer  model  used  is  con¬ 
sidered  adequate  for  the  purpose  of  this  investigation. 

From  the  comparison  of  torque  and  leakage,  it  is  known  that 
using  a  3  0.1  un  and  F  =  0.15  provides  a  reasonable  prediction  of 
performance.  Thus,  a  set  of  load  versus  film  thickness  curves 
using  a  -  0.1  nn  and  F  *  0.15  should  provide  some  prediction 
about  puffing  of  the  seal.  These  curves  are  shown  in  Figure 
4-26.  The  progression  of  curves  shows  that  for  B  =  1.0  seal 
there  are  three  roots  for  the  T*  3  445  *K  curve.  However,  they 
are  close  together  such  that  puffing  would  not  occur.  This 
agrees  with  the  experimental  result. 

As  for  B  3  0.75  seal,  considering  the  curve  of  T*,  =  472*K 
( 199* C) ,  although  there  are  three  roots,  the  right-hand  root  does 
not  move  very  far  to  the  right  before  disappearing.  This  means 
that  noticeable  puffing  may  not  occur  for  B  =  0.75  seal  either. 
However,  puffing  has  been  observed  at  T„  =  464*K  (191*C)  during 
Test  73.  One  possible  explanation  for  this  is  that  the  balance 
ratio  of  the  pressurized  seal  ring  might  change  due  to  mechanical 
and  thermal  distortion  in  operation.  If  it  changed  to  as  little 
as  B  =*  0.74,  then  the  horizontal  line  of  applied  load  would  be 
lower  and  the  right-hand  root  would  move  far  to  the  right.  This 
would  cause  puffing. 

Therefore,  based  on  the  comparison  above,  the  theoretical 
model  cannot  provide  an  accurate  prediction  for  puffing.  How¬ 
ever,  the  model  does  have  some  use  in  that  it  indicates  the 
approximate  temperature  at  which  puffing  would  be  expected  to 
occur. 

Both  the  torque  versus  temperature  and  aT  versus  temperature 
data  display  a  wide  scatter  in  spite  of  being  averaged  over  many 
hours  of  testing.  As  shown  in  Chapter  2,  numerous  attempts  were 
made  to  get  repeatable  results  at  a  given  temperature.  This 
scatter  is  similar  to  scatter  found  in  all  previous  tests  for 
para1 lei  face  seals.  It  is  believed  that  certain  random 
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Figure  4-26.  Load  Function. 


processes  are  occurring  to  cause  this.  The  deterministic  model 
used  to  predict  performance  is  simply  not  adequate.  It  is 
thought  that  there  is  an  interaction  between  wear  (global  and 
local),  thermal  deformation  (global  and  local),  friction  coeffi¬ 
cient  (local)  and  fluid  flow  (local)  which,  in  fact,  may  allow  a 
seal  to  operate  at  different  average  conditions  from  time  to  time. 
Clearly,  more  investigation  of  the  basic  friction  and  wear  proces¬ 
ses  is  needed  to  explain  this  type  of  behavior. 

Offset-Land  Seal 

Using  the  methods  outlined  in  the  previous  chapter,  perform¬ 
ance  was  calculated  for  the  experimental  offset-land  seal  de¬ 
scribed  in  Chapter  2.  The  computed  results  are  shown  compared  to 
experimental  results  in  Figure  4-27.  Only  two  experimental 
points  were  plotted  because  in  only  two  of  the  tests  could  the 
final  offset  be  determined.  It  is  no  surprise  that  one  of  the 
torque  readings  is  poorly  predicted  since  the  torque  values  are  at 
the  limit  of  accurate  resolution  of  the  torque  measuring  device. 
Considering  leakage,  the  theory  underpredicts  at  the  lower  offset 
case.  This  probably  occurs  because  the  roughness  of  the  surface 
of  the  land  was  larger  than  assumed  in  the  calculations. 

In  summary,  some  agreement  between  theory  and  experiment  is 
observed.  Clearly,  many  more  experimental  results  will  be  needed 
before  any  final  judgement  can  be  made  concerning  the  validity  of 
the  theory. 

Speed  Effects 

In  the  previous  report  [8],  computations  of  the  speed  ef¬ 
fects  were  made  and  compared  to  experimental  results.  These  cal¬ 
culations  have  now  been  remade  using  the  new  algorithm  described 
previously  and  compared  to  all  of  the  speed  effect  data  now  avail¬ 
able.  Tables  4-4  and  4-5  show  the  computed  values  assuming  0% 
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wear  and  100%  wear,  respectively.  The  0%  wear  case  represents 
expected  performance  if  little  wear  occurred  and  should  predict 
performance  at  start  up  or  under  conditions  where  the  worn  shape 
of  the  face  was  not  fully  developed.  The  100%  wear  case  repre¬ 
sents  expected  performance  for  a  well  run  in  case. 

Figures  4-28  and  4-29  compare  the  theoretical  and  experi¬ 
mental  results.  All  of  the  speed  effect  experimental  results  for 
n  *  3  are  plotted.  Considering  torque,  agreement  with  the  100% 
wear  case  at  and  above  900  rpm  is  satisfactory.  At  0  rpm,  exper¬ 
imental  torque  is  much  higher  than  predicted  by  even  the  0%  wear 
case.  Furthermore,  experimental  results  consistently  show  a  tor¬ 
que  minimun  around  450  rpm,  whereas  theory  shows  the  minimum  at 
* 

just  above  900  rpm.  At  these  lower  speeds  the  theory  is  clearly 
missing  some  important  physical  behavior.  One  good  possibility 
is  that  the  assumed  mechanical  coefficient  of  friction  is  not 
constant  but  may  drop  off  quickly  with  increasing  speed  due  to 
microasperity  lubrication  effects.  That  is,  at  the  lowest  speeds 
tested  (1.4  rpm)  microasperity  lubrication  does  not  function  and 
the  mechanical  coefficient  of  friction  is  much  higher  than  as¬ 
sumed.  Such  behavior  could  well  explain  the  discrepancy  in  the 

/ 

torque  prediction,  and  this  question  will  probably  become  re¬ 
solved  as  better  data  for  mechanical  friction  coefficients  in 
water  become  available  using  the  friction  test  apparatus  de¬ 
scribed  in  Chapter  2. 

Considering  leakage  in  4-29,  disagreement  at  the  lower 
speeds  is  pronounced.  It  is  probable  that  at  the  lower  speeds 
100%  wear  has  not  been  achieved  and  the  lower  leakage  values 
might  correspond  better  to  the  0%  wear  case  which  they  do.  It  is 
also  possible  that  there  are  additional  deflection  components  of 
the  carbon  that  have  not  been  accounted  for.  Finally,  it  is 
likely  that  the  surface  roughness  used  in  the  calculation  does 
not  accurately  characterize  the  actual  roughness  as  discussed 
previously. 
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In  any  case  the  wavy  seal  theory  appears  to  be  inadequate  at 
the  lower  speeds  even  though  Figures  4-1,  4-2,  4-28  and  4-29  show 
reasonable  agreement  at  higher  speeds.  It  is  thought  that  some 
insight  into  these  problems  may  arise  from  data  from  a  new  nine- 
wave  optimum  design  which  has  negligible  deflection  and  from  the 
friction  and  wear  test  apparatus.  As  new  information  becomes 
available,  the  wavy  seal  model  will  be  modified  accordingly. 

Nine-Wave  Seal 

As  described  in  Chapter  2  tests  were  run  using  the  original 
three-wave  seal  as  a  nine-wave  seal.  Resu.ts  using  the  wavy  seal 
model  have  been  computed  and  are  shown  in  Table  4-6.  Friction 
torque  is  seen  to  be  consistently  high  at  all  speeds  as  it  was  in 
the  tests  and  leakage  is  negligibly  small  as  in  the  tests.  The 
calculations  show  that  this  particular  nine-wave  seal  performs 
very  poorly  as  it,  in  fact,  did. 

It  is  interesting  to  note  that  using  the  original  seal  model 
algorithm  a  good  performance  was  predicted.  This  occurred  be¬ 
cause  the  original  algorithm  could  not  create  the  lubricant  star¬ 
vation  that,  in  fact,  occurs  as  shown  by  the  Elrod  algorithm  seal 
model.  In  this  case,  lubrication  starvation  is  caused  by  the 
fact  that  the  effectiveseal ing  radius  is  predicted  to  be  near 
the  outside  diameter  and  not  close  to  the  inside  diameter  as  it 
must  be  to  approach  optimum  performance.  Because  of  this  an 
insufficient  amount  of  fluid  is  present  in  the  highly  converging- 
diverging  region  of  the  seal  to  provide  a  significant  hydrodyna¬ 
mic  effect.  In  the  original  algorithm,  this  did  not  occur  be¬ 
cause  negative  pressures  were  allowed  to  pull  fluid  in  from  the 
low  pressure  side. 

In  summary,  the  seal  model  predicts  the  poor  measured  per¬ 
formance  of  this  particular  nine-wave  seal.  This  example  also 
serves  to  illustrate  how  the  original  seal  model  based  on  approx¬ 
imate  cavity  boundary  conditions  could  lead  to  very  poor  design 
because  of  its  inaccuracy. 


150 


Chapter  5 

Sumary,  Conclusions,  and  Recommendations 


Experimental  Results 

1)  The  waviness  drive  provides  a  means  of  producing  a  mov¬ 
ing  wave  using  no  moving  parts  internal  to  the  seal.  The  device 
designed  performed  with  no  problems  for  a  500-hour  test.  Seal 
performance  with  time  does  not  deteriorate  using  this  device  as 
compared  to  sliding  0-rings.  Minor  modifications  and  improve¬ 
ments  of  the  device  to  smooth  out  the  pressure  versus  time  func¬ 
tions  are  anticipated. 

2)  For  the  purpose  of  obtaining  friction  and  wear  data  in 
water  but  independently  of  hydrostatic  pressure,  a  friction  and 
wear  test  apparatus  has  been  designed.  It  is  thought  that  data 
obtained  from  this  device  will  explain  some  of  the  discrepancies 
between  theory  and  experiment  found  using  the  simple  wear  and 
friction  relationships  built  into  the  present  seal  models. 

3)  Experimental  results  from  a  series  of  tests  at  different 
radial  tapers  show  that  thermally  caused  radial  taper  has  a 
strong  influence  on  initial  seal  performance.  A  large  initial 
divergent  taper  must  be  lapped  onto  the  seal  faces  such  that  ini¬ 
tial  operation  will  have  the  faces  parallel.  Results  also  show 
that  initially  convergent  faces  result  in  very  low  friction  but 
some  wear  at  the  inside  corners  occurs.  This  will  eventually 
lead  to  parallel  face  operation.  However,  the  amount  of  time 
required  to  establish  parallel  face  operation  after  starting  at  a 
high  convergent  value  is  much  larger  than  the  few  hundred  hours 
of  test  used  herein.  Initial  seal  wear  can  definitely  be  reduced 
by  starting  with  a  convergent  taper.  The  disadvantage  of  doing 
this  is  that  long-term  performance  cannot  be  anticipated  by 
short-term  testing.  It  is  thought,  in  fact,  that  long-term  per¬ 
formance  can  best  be  simulated  by  lapping  a  divergent  radial 
taper  equivalent  to  the  thermal  taper  onto  the  seal.  It  was 
found  that  the  behavior  of  a  balanced  seal  operating  as  described 
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was  much  different  from  seals  where  no  significant  initial  radial 
taper  was  provided.  Thus,  testing  seals  with  no  initial  taper 
can  be  misleading  because  of  the  thermal  rotation  effect. 

4)  Results  from  the  high  temperature  tests  show  that  as  the 
seal  operating  temperature  goes  above  the  saturation  temperature 
of  the  liquid  at  discharge  pressure,  seal  friction  is  reduced. 
This  trend  continues  until  just  before  saturation  temperature  at 
the  seal  inlet  pressure  is  reached.  At  this  point,  friction 
greatly  increases  with  increasing  temperature.  There  is  a  great 
amount  of  scatter  in  the  data  which  cannot  be  attributed  to 
experimental  error.  It  appears  possible  that  multiple  states  of 
equilibrium  can  be  set  up  for  a  high  temperature  parallel  face 
seal  just  as  for  a  normal  temperature  parallel  face  seal.  The 
question  of  why  the  high  temperature  seals  formed  a  converging 
final  taper  instead  of  the  normal  diverging  taper  remains 
unanswered . 

5)  The  radial  taper  tests,  the  high  temperature  tests,  and 
parallel  face  tests  show  that  for  parallel  faces  the  use  of  the 
Patir-Cheng  [16]  surface  roughness  model  to  calculate  leakage 
overestimates  the  leakage  flow  by  at  least  an  order  of  magni¬ 
tude.  Investigation  of  the  worn  carbon  surface  shows  that  the 
surface  consists  of  mostly  high  regions  having  a  small  roughness 
and  deep  craters  which  are  generally  not  interconnected.  Thus, 
while  the  craters  add  to  the  roughness  value  of  the  surface,  they 
do  not  promote  flow.  It  is  clear  that  the  surfaces  are  not 
Gaussian  as  assumed  for  the  purpose  of  using  the  model  and  have 
some  very  special  function.  More  work  needs  to  be  done  to  char¬ 
acterize  the  flow  between  the  rough  surfaces  of  seals.  It  is 
thought  that  a  model  similar  to  that  of  Patir  and  Cheng  could  be 
developed  based  on  actual  surface  characteristics.  A  better 
rough  surface  flow  model  is  needed  to  use  in  the  various  other 
seal  models  to  improve  their  predictive  capabilities. 


154 


6)  The  results  from  two  500-hour  test  show  that  the  wavi- 
ness  concept  reduces  wear  by  a  minimun  of  a  factor  of  200  com¬ 
pared  to  a  flat-faced  test.  Wear  is  so  small  that  it  is  of  the 
same  order  of  magnitude  as  the  carbon  nose  distortions  and  some¬ 
times  negative  wear  values  are  measured.  Improved  wear  measuring 
techniques  are  being  incorporated  into  future  seal  tests. 

7)  The  offset  land  seal  concept  works  in  principle.  How¬ 
ever,  it  is  not  known  if  the  protrusions  will,  in  fact,  wear  away 
over  a  long  period  of  operation.  It  was  not  possible  to  start  a 
seal  in  the  offset  low  friction  mode  at  low  values  of  offset.  At 
high  values  of  offset  the  seal  works  good  but  leakage  is  high. 
Thus,  at  this  point  it  is  not  knovwi  particularly  if  the  offset 
concept  can  be  made  to  work  to  provide  a  low  leakage,  low  fric¬ 
tion  seal  as  originally  intended.  More  experimental  work  parti¬ 
cularly  at  low  offsets  needs  to  be  performed. 

8)  Additional  data  was  obtained  at  various  speeds  for  both 
the  wavy  seal  and  the  parallel  face  seal.  At  very  low  speeds 
friction  is  quite  large  in  both  cases.  There  is  a  rapid  decrease 
in  friction  with  increasing  speed  for  both  the  wavy  and  parallel 
face  seals.  The  wavy  seal  friction  is  reduced  much  further  by 
speed  increases  than  the  parallel  face  seal.  There  is  strong 
evidence  of  micro asperity  lubrication  in  both  cases,  that  is, 
friction  decreases  with  increasing  speed  when  no  macro  scale 
hydrodynamic  lift  is  present.  It  is  expected  that  the  use  of  the 
wear  test  apparatus  will  help  to  resolve  this  question. 

9)  Tests  using  the  existing  seal  as  a  nine-wave  seal  show 
poor  performance.  Friction  is  high.  Later  theory  shows  this  to 
be  expected  and  indicates  that  a  wavy  seal  cannot  be  casually 
designed  and  expect  good  performance.  Optimum  wavy  seal  design 
principles  must  be  adhered  to. 

Theoretical  Results 

1)  The  principles  governing  the  waviness  generator  are 


quite  simple  and  similar  to  those  for  a  poly  phase  electric 
motor.  To  this  author's  knowledge,  this  principle  has  not  been 
previously  applied  to  a  mechanical  device.  There  may  be  other 
applications  for  this  principle. 

2)  The  seal  waviness  model  has  been  reworked  using  the 
Elrod  algorithm  to  predict  cavitation.  Although  compared  to  the 
previous  model,  little  difference  in  prediction  was  observed  for 
the  three-wave  seals,  a  great  difference  was  noted  for  the  non¬ 
optimum  nine-wave  seal  above.  The  previous  model  overestimates 
load  support.  The  Elrod  algorithm  based  seal  model  is  therefore 
more  accurate  and  should  be  used  hereafter.  A  small  increase  in 
computer  time  is  needed  for  the  new  seal  model,  but  no  serious 
limitation  is  imposed. 

3)  A  paper  was  written  on  optimun  wavy  seal  design  [14]. 

It  has  been  shown  that  for  seals  operating  under  variable  speeds 
and  pressures,  a  relatively  large  number  of  waves  must  be  used  so 
that  the  relative  tilt  of  the  faces  does  not  change  under  vari¬ 
able  conditions.  In  the  three-wave  seal,  results  show  that 
deflections  change  with  pressure.  This  will  lead  to  accelerated 
wear.  It  has  been  determined  that  about  nine  waves  are  required 
for  sufficient  stiffness  for  seals  of  the  size  herein.  Further 
refinement  in  determining  the  optimun  wavy  seal  configuration 
beyond  the  discussion  in  the  previous  report  [8]  has  been  pro¬ 
vided.  A  new  seal  based  on  these  principles  is  being  designed. 

4)  The  theory  supporting  the  offset  land  concept  is  devel¬ 
oped.  Theory  shows  that  an  offset  land  seal  with  very  low  fric¬ 
tion,  wear  and  leakage  can  be  designed. 

Comparison  of  Theoretical  and  Experimental  Results 

1)  Comparison  of  theory  and  experiment  for  three-wave  seals 
at  6.9  MPa  gas  pressure  is  good. 

2)  The  radial  taper  model  provides  a  good  prediction  for 
the  short-term  effects  of  varying  the  radial  taper.  The  effect 
of  thermal  rotation  is  properly  accounted  for.  The  long-term 
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performance  is  poorly  predicted.  It  is  thought  that  the  wear 
model  used  is  inadequate  and  it  is  expected  that  wear  test 
results  obtained  on  the  friction  and  wear  test  apparatus  will 
lead  to  improved  predictions. 

3)  The  high  temperature  seal  model  predicts  the  trend  and 
transition  point  of  the  high  temperature  tests.  The  model  cannot 
predict  or  explain  the  scatter  of  the  data.  It  is  believed  that 
certain  random  processes  may  account  for  this  scatter.  A  sto¬ 
chastic  seal  model  is  needed  to  better  understand  this  behavior. 

4)  Some  agreement  between  theory  and  experiment  is  observed 
for  the  offset  land  studies.  More  experimental  data  is  needed  to 
refine  this  comparison. 

5)  While  predicting  the  general  trend,  the  wavy  seal  model 
does  a  poor  job  on  predicting  wavy  seal  performance  as  a  function 
of  speed.  It  is  thought  that  refined  friction  measurements  will 
be  helpful.  However,  it  appears  that  other  factors  such  as  seal 
deflection  may  also  be  involved. 

6)  The  refined  wavy  seal  model  does  predict  the  high  fric¬ 
tion  and  low  leakage  observed  for  the  nonoptimum  nine-wave  seal. 

General 

The  experimental  and  theoretical  foundations  for  the  moving 
wave  seal  have  been  well  established.  It  has  been  shown  experi¬ 
mentally  under  somewhat  ideal  conditions  that  wear  is  insignifi¬ 
cant.  It  has  been  shown  that  a  waviness  device  and  seal  can  be 
operated  for  a  500  hour  test  with  no  mechanical  problem  and  no 
deterioration  of  purpose.  It  has  been  shown  that  the  wavy  seal 
model  does  a  reasonable  job  of  predicting  performance  although 
certain  improvements  are  needed.  In  particular,  a  better  under¬ 
standing  of  the  effect  of  actual  surface  roughness  on  leakage, 
wear  mechanisms,  and  the  change  of  friction  coefficient  under 
various  operating  conditions  are  needed.  Experimental  programs 
have  been  outlined  to  provide  this  needed  information.  Variables 
which  govern  optimun  wavy  seal  performance  are  well  understood. 
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Test  #53,  T  =  37.8°C,  pH  n  =  1.72  MPa,  1800  RPM 
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Test  #54,  T*  =  93.3°C,  p„  n  =  1.72  MPa,  1800  RPH 
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Test  158,  T.,,  =  176. 7°C,  pH  ft  =  1.72  MPa,  1800  RPH 
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Test  §62 ,  pu  n  =  l-72  MPa»  1800  RpM 
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Test  #64.  1800  RPM,  pH  n  =  3.45  MPA,  =  1013  Mm/m  ,  B  =  0.75 
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Test  *72,  T„  =  37.8°C,  pu  n  =  1.72  MPa,  1800  RPM,  B  =  0.75 
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Test  #77.  1800  RPM,  PH  n  =  3.45  MPa,  P„  =  5.52  MPa 
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Test  #94.  900  RPM,  Pu  =  3.45  MPa,  P  =  0.0 
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Test  #98.  1800  RPM,  P„  n  =  3-45  MPa,  =  1200  wm/m,  B  =  0.75 


TEST#  99 
10/23/80 


TEST#  100 


Test  #100.  1800  RPM,  P„  n  »  3.45  MPa,  =  643  pm/m,  B  =  0.75 


TEST#  101 
10/28/80 


TIME  (h) 

Test  #101.  1800  RPM,  P„  «  *  3.45  MPa,  =  660  ym/m,  B  =  0.75  ,  50%  wear 


TIME  (h) 


(UMu/guiO)  33BHB31 

Ln^r^rfnnruru  —  — h 


m  (S 


LU 

z: 


CO 


CO 


(3.)  3BniBB3dW3i  “IB3S 


246 


Test  #103.  1.4  RPM,  P„  n  =  3.45  MPa,  P  =  6.9  MPa 


(o,)  3aruua3dW3i  ito 

247 


TEST#  105 
11/13/80 
Temperature 


TIME  (h) 

Test  #105.  1800  RPM,  P„  n  =  3.45  MPa,  B  =  0.75  (Land  Seal) 


TEST#  106 
11/17/80 


Test  #106.  1800  RPM,  PH  n  =  3.45  MPa,  B  =  0.75  (Land  Seal) 


TIME  (h) 


TEST*  108 
12/3/80 


TEST#  109 
1/8/81 


Efrect  of  waviness  pressure  on  torque 


TIME  Ch) 

Test  #110.  1800  RPM,  PH  =  3.45  MPa,  P  =  6.9  MPa,  n  =  3  (with  new  waviness  generator) 


TEST*  110 
1/25/81 


CD 

in 


(Do)  3aruya3dW3i  ib3S 
256 


TIME  (h) 
(continued) 


TIME  <h) 

(continued) 


Distribution  List 


Recipient 

Office  of  Naval  Research 

800  N.  Quincy  Street 

Arlington,  Virginia  22217 

Attn:  M.  Keith  Ellingsworth,  Code  473 

Defense  Documentation  Center 
Building  5 
Cameron  Station 
Alexandria,  Virginia  22314 

Naval  Research  Laboratory 
4555  Overlook  Avenue 
Washington,  DC  20390 
Attn:  Technical  Information  Division 
Code  2627 

Dr.  Ravner,  Code  6170 

U.S.  Naval  Postgraduate  School 

Monterey,  California  93940 

Attn:  Dept,  of  Mechanical  Engineering 

U.S.  Naval  Academy 

Annapolis,  Maryland  21402 

Attn:  Dept,  of  Mechanical  Engineering 

Naval  Air  Systems  Command 
Jefferson  Plaza 
Washington,  DC  20360 
Attn:  B.  Poppert,  Code  240E 

Naval  Sea  Systems  Command 
Crystal  City,  National  Center  #3 
Washington,  DC  20360 
Attn:  Frank  Ventriglio,  Code  0SR14 

Naval  Ships  R&D  Center 
Annapolis,  Maryland  21402 
Attn:  Friction  and  Wear  Branch 
J.  F.  Dray 

Naval  Air  Engineering  Center 
Lakehurst,  New  Jersey  08733 
Attn:  Mr.  P.  Senholzi 


Number  of  Copies 


(3) 


(12) 


(1) 

0) 

(1) 


(1) 


(1) 


(1) 

O) 


Naval  Air  Propulsion  Test  Center 
Trenton,  New  Jersey  08628 
Attn:  Mr.  R.  Valori 


259 


flWBB*  PAK  KdkMMJOT 


Distribution  List  (continued) 


Recipient 

Naval  Air  Development  Center 
Warminster,  Pennsylvania  18974 
Attn:  Mr.  A.  Conte 

National  Science  Foundation 
1800  G  Street,  NW 
Washington,  DC  20550 
Attn:  Dr.  C.  J.  Asti  11 

National  Bureau  of  Standards 
Washington,  DC  20234 
Attn:  Dr.  W.  Ruff 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  Ohio  44135 
Attn:  R.  L.  Johnson 

Air  Force  Office  of  Scientific  Research 
Washington,  DC  20333 

Attn:  Directorate  of  Engineering  Sciences 

Air  Force  Aeropropul sion  Laboratory 
Wright-Patterson  Air  Force  Base, 

Ohio  45433 

Attn:  AFAPL/P0D-1 ,  Dick  Quigley,  Jr. 

Army  Research  Office 
Durham,  North  Carolina  27706 
Attn:  Dr.  £.  A.  Saibei 

Office  of  Naval  Research  Branch  Office 
1030  East  Green  Street 
Pasadena,  California  91106 

Assistant  Chief  for  Technology 
Office  of  Naval  Research,  Code  200 
Arlington,  Virginia 

Prof.  H.  S.  Cheng 

Department  of  Mechanical  Engineering 
Northwestern  University 
Evanston,  Illinois 

Crane  Packing  Company 
6400  Oakton  Street 
Morton  Grove,  Illinois  60053 
Attn:  Art  Zobens 


Number  of 
(1) 

(1) 


(1) 

(1) 


(1) 

(1) 


(1) 

(1) 

(1) 

(1) 


(1) 


260 


Distribution  List  -  continued 


Recipient 

Sealol,  Inc. 

Box  2158 

Providence,  Rhode  Island  02905 
Attn:  H.  F.  Greiner 

Pure  Carbon  Company 

St.  Marys,  Pennsylvania  15857 

Attn:  R.  R.  Paxton 

Franklin  Research  Institute 
20th  and  Race  Streets 
Philadelphia,  Pennsylvania  19103 
Attn:  Harry  C.  Rippel 

Naval  Sea  Systems  Command 
Crystal  City,  National  Center  #3 
Washington,  DC  20360 
Attn:  Code  524,  Dick  Graham 

Naval  Ships  R&D  Center 
Annapolis,  Maryland  21402 
Attn:  A1  Harbaugh 

Prof.  Francis  E.  Kennedy,  Jr. 

Thayer  School  of  Engineering 
Dartmouth  College 
Hanover,  New  Hampshire  03755 

Mr.  Cl  if  Mussen 
PMS  396-223 

Trident  Submarine  Ship  Acquisition  Project 
Naval  Sea  Systems  Command 
Department  of  Navy 
Washington,  DC  20362 


261 


